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1  Introduction 


Problem  Description 

Outlet  works  consist  of  a  combination  of  structures  designed  to  control  the 
release  of  water  from  the  reservoir  as  required  for  project  purposes  or  operation. 
Outlet  works  structures  are  typically  found  in  embankment  dams,  although  in 
some  cases  concrete  dams  may  exhibit  outlet  works  detached  from  the  dam.  Out¬ 
let  works  are  used  to  pass  diversion  flows  during  constmction,  regulate  flood 
flows,  support  water  supply  and  power  generation  operations,  aid  in  emptying  the 
reservoir  in  an  emergency  condition,  and  permit  reservoir  lowering  for  inspec¬ 
tions  and  special  repairs.  In  many  dams,  outlet  works  provide  the  only  means  of 
lowering  the  reservoir  and  thereby  providing  a  means  of  minimizing  the  risk  of 
catastrophic  failure  in  the  event  of  an  earthquake.  Their  proper  functioning 
immediately  after  an  earthquake  is  very  important,  since  the  controlled  release  of 
the  reservoir  could  help  to  prevent  additional  propagation  of  damage  by  reducing 
the  hydrostatic  loads. 

Typically,  outlet  works  consist  primarily  of  an  intake  structure  and  a  conduit 
or  tunnel,  in  addition  to  other  features  necessary  for  operation,  such  as  control 
gates,  approach  and  discharge  channels,  and  energy  dissipaters  (Headquarters, 
U.S.  Army  Corps  of  Engineers  2003).  The  intake  structure  may  exhibit  multiple 
openings  with  die  corresponding  valves  or  gates,  to  allow  for  more  effective  con¬ 
trol  of  the  flow  of  water.  The  control  mechanisms  are  typically  located  within  the 
intake  tower,  although  configurations  with  control  gates  within  the  conduit  or  at 
the  downstream  end  of  the  conduit  are  also  possible.  Intake  towers  can  be  sepa¬ 
rated  into  two  broad  categories:  free-standing  (vertical  structure,  typically  with  a 
service  bridge  at  the  top)  and  inclined  (stracture  supported  against  the  abutment). 
This  research  work  will  concentrate  on  free-standing  intake  towers,  since  they 
constitute  the  most  common  structural  configuration. 

The  dynamic  response  of  an  intake  tower  during  an  earthquake  may  present 
quite  complex  characteristics  due  to  many  factors.  The  water,  inside  and  outside 
of  the  tower,  plays  an  important  role  in  the  modification  of  the  response  of  the 
structure.  In  addition,  the  complexity  of  the  analysis  is  increased  by  important 
soil-stracture  interaction  effects.  In  addition,  although  this  effect  is  generally 
neglected,  the  intake  structure  is  actually  subjected  to  multiple-support  excita¬ 
tions.  Although  a  free-standing  intake  tower  is  basically  a  hollow  column,  it 
exhibits  an  inspection  bridge  that  connects  the  top  of  the  tower  to  the  dam  or 
abutment.  Therefore,  the  intake-outlet  tower  and  bridge  system  has  two 
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well-separated  supporte.  Different  local  soil  conditions,  teveling  wave  effects, 
phase  difference  in  the  seismic  motions,  and  other  fectors  can  cause  a  spatial 
variability  of  the  ground  motion  at  the  supports.  These  multiple-support  excita¬ 
tion  effects  can  be  important  in  Ihe  calculation  of  the  seismic  response  of  die 
system.  Therefore,  and  because  of  the  lack  of  previous  investigations  dealing 
with  these  effects,  there  is  a  need  to  study  in  more  detail  their  influence  on  the 
djmamic  response  of  the  intake  tower. 

The  main  objective  of  die  research  described  herein  was  to  calculate  the 
dynamic  response  of  the  intake-outlet  tower  under  earthquake  excitations,  giving 
proper  consideration  to  multiple-support  effects  and  soil-stracture  interaction 
phenomena.  Several  alternative  methods  were  considered  for  the  analyses  in 
Older  to  compare  the  results  obtained  with  a  simplified  model  of  the  tower-bridge 
system  with  those  generated  using  more  refined  finite  element  models.  It  is 
advantageous  to  show  that  the  resulte  from  a  simplified  model  are  comparable  to 
those  from  the  more  accurate  finite  element  models,  since  practicing  engineers 
will  in  all  likelihood  make  use  of  the  simplified  model. 


Justification 

Several  natural  hazards  affect  both  human  hves  and  properties.  Earthquake 
ground  motions  are  one  of  the  most  destractive  hazards,  and  this  threat  is 
increased  by  its  random  nature.  Earthquakes  affect  all  human-made  stractures 
founded  on  ground  or  buried.  Ihe  feilure  of  certain  structures  also  compromises 
the  lives  of  many  people  living  in  the  surrounding  areas.  Dams  are  among  these 
critical  structures;  thus,  they  have  been  and  remain  the  object  of  attention  of  the 
stractural  and  geotechnical  engineering  community.  Intake-outlet  towers  could 
play  an  important  role  in  preventing  the  catastrophic  failure  of  a  dam  after  an 
earthquake,  by  providing  a  means  of  controlling  the  water  level  and  reducing  the 
corresponding  hydrostatic  pressure.  Ihe  stractural  safety  of  these  towers  and 
their  fimctional  capability^  in  the  event  of  a  major  earthquake  are  crucial. 

As  pointed  out  by  Dove  (1996),  many  intake  towers  underfire  responsibility' 
of  the  U.S.  Army  Ckrrps  of  Engineers  (USAGE)  are  located  in  areas  associated 
with  high  seismic  activity.  Therefore,  fire  analysis  of  the  stractural  safety'  and 
functional  c^ability  of  these  towers  in  the  event  of  a  major  earthquake  reaches 
paramount  importance.  Observations  during  past  earthquakes  demonstrated  the 
importance  of  the  effects  of  multiple  excitations  in  some  extended  stractures, 
such  as  bridges,  pipelines,  and  dams  (Der  Kiureghian  and  Neuenhofer  1991).  In 
studies  that  are  currently  performed  for  the  analysis,  design,  and  verification  of 
intake  towers,  these  mirltiple  excitation  effects  are  not  considered.  This  feet 
increases  the  need  to  develop  better  procedures  for  the  accurate  prediction  of  the 
dynamic  response  of  these  structures.  This  report  intends  to  address  this  concern 
by  developing  a  relatively  simple  model  to  perform  seismic  response  analysis  of 
free-standing  intake  towers  considering  the  spatial  variability  of  the  ground 
motion  and  soil-stracture  interaction  effects.  In  aidition,  this  report  will  evaluate 
the  importance  of  considering  these  tw'o  effects  in  the  analysis. 
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Contribution  to  the  State  of  the  Art 

To  the  authors’  knowledge,  all  previous  studies  dealing  with  free-standing 
intake  towers  did  not  consider  the  effects  of  attached  structures  such  as  the  access 
bridge  (Figure  1),  and  the  models  used  in  the  analyses  were  limited  to  the  towers 
themselves.  In  general  engineering  practice,  a  typical  “stick  model”  is  used  to 
study  these  towers  when  the  geometry  of  the  tower  is  relatively  simple  and  the 
bridge  has  a  light  weight.  However,  it  is  reasonable  to  expect  that  the  bridge  can 
have  an  effect  on  the  response  of  the  tower.  Whether  the  effect  of  the  bridge  is 
such  that  it  reduces  or  magnifies  the  seismic  response  is  not  clear.  Thus,  a  model 
of  the  bridge  will  be  included  to  calculate  the  seismic  response  in  order  to  assess 
its  importance. 


Figure  1 .  Model  of  the  intake-outlet  tower  plus  the 
access  bridge 


The  inclusion  of  a  model  of  a  bridge  can  be  easily  accomplished  either  by 
modeling  the  bridge  as  a  flexural  beam  that  is  rigid  in  the  axial  direction  or  by 
using  frame  elements  including  axial  deformation.  Appropriate  definition  of  the 
ground  acceleration  is  more  complex.  It  is  reasonable  to  expect  that  the  absolute 
acceleration  at  the  base  of  the  tower  may  not  be  the  same  as  that  experienced  by 
the  tower  at  the  elevation  where  the  bridge  connection  is  made.  Therefore,  the 
effects  of  traveling  waves,  lack  of  coherence,  and  other  phenomena  associated 
with  multiple-support  excitations  must  be  considered. 

One  could  argue  that  these  effects  may  not  be  important  since  the  distances 
between  the  supports  of  the  tower  and  bridge  are  not  substantial.  This  may  be 
true  in  some  cases,  but  the  fact  that  the  supports  are  at  different  elevations  can 
heighten  the  dissimilar  character  of  the  ground  acceleration.  The  ground  accele¬ 
ration  could  suffer  amplification  or  attenuation  depending  on  the  soil  properties. 
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the  height  of  the  dam,  and  the  frequency  content  of  the  excitation.  The  previous 
studies  dealing  with  multiple-support  excitations  considered  that  the  supports  of 
the  extended  structure  were  well  separated  in  space  but  at  the  same  elevation. 
This  is  the  case  of  the  highway  bridges  and  pipeline  systems.  Therefore,  one  of 
the  areas  that  will  be  studied  is  the  case  in  which  the  supports  are  at  different 
levels.  In  addition,  other  analyses  will  be  performed  considering  the  soil-structure 
interaction  effects  to  obtain  a  better  insight  into  the  dynamic  response  of  this  type 
of  structure. 

The  first  step  undertaken  to  consider  the  multiple-support  excitation  will  be 
to  prepare  a  finite  element  model  of  the  topography  of  the  earthen  dam.  This 
model  will  be  used  to  calculate  the  response  at  the  location  of  tiie  bridge.  The 
records  of  the  absolute  acceleration  at  the  latter  point  and  at  the  base  are  used  as 
input  for  a  combined  model  of  the  bridge  and  tower.  A  computer  program  will  be 
developed  to  account  for  the  effects  of  the  multiple-support  excitations  and  the 
soil-structure  interaction  in  the  tower  and  access  structure. 


Report  Organization 

Chapter  1  contains  a  general  introduction.  The  justification  of  the  research 
and  the  problem  description  are  briefly  discussed.  The  chapter  includes  a  discus¬ 
sion  of  the  specific  contribution  to  the  state  of  the  art. 

Chapter  2  contains  a  brief  literature  review  related  to  the  topics  of  the 
research.  The  first  section  of  this  chapter  presents  information  regarding  studies 
of  intake-outlet  towers.  The  seismic  response  of  these  towers  is  discussed  in  the 
second  section.  This  section  explains  the  factors  that  affect  the  response  of  this 
type  of  structure,  such  as  the  fluid-structure  and  soil-structure  interaction  and 
multiple-support  excitation  effects.  The  third  section  is  devoted  to  the  general 
concepts  and  a  review  of  the  most  relevant  works  on  the  subject  of  multiple- 
support  excitation  effects.  Also  included  is  a  brief  review  of  several  investiga¬ 
tions  related  to  soil-structure  interaction  effects  in  structures  such  as  intake-outlet 
towers. 

Chapter  3  describes  the  general  methodology  used  to  accomplish  the  goals  of 
this  report.  The  chapter  starts  with  an  introduction  of  the  assumptions  and  consid¬ 
erations  taken  into  account  during  the  research.  There  is  a  brief  description  of  the 
finite  element  models  and  computer  programs  developed  for  this  research.  There 
is  also  a  description  of  the  procedure  followed  to  select  the  specific  tower  used  in 
the  study.  Drawings  showing  the  final  intake-outlet  tower  selected  for  the  case 
study  are  presented  in  this  section. 

Chapter  4  is  devoted  to  the  description  of  the  finite  element  models  devel¬ 
oped  during  this  research  and  previously  mentioned  in  the  general  methodology 
(Chapter  3).  The  chapter  includes  descriptions  of  the  SAP2000  and  QUAD4M 
computer  programs,  which  are  the  principal  tools  for  the  finite  element  analyses 
used  in  this  research.  The  rest  of  the  chapter  contains  detailed  information 
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regarding  the  finite  element  models  developed  during  the  study,  such  as  dimen¬ 
sions,  boundary  conditions,  and  input  data. 

Chapter  5  contains  the  description  of  the  computer  program  developed  in  the 
MATLAB  computational  environment,  which  is  the  main  tool  used  in  this 
research.  The  formulations  upon  which  this  computer  program  is  based  are  dis¬ 
cussed  in  several  sections.  First,  the  multiple-support  excitation  formulation  is 
presented  in  detail.  The  simplified  approach  developed  by  Goyal  and  Chopra 
(1989a)  for  the  calculation  of  the  added  hydrodynamic  masses  along  the  height  of 
the  tower  is  also  included.  The  procedure  to  account  for  the  soil-structure  inter¬ 
action  effects  in  the  towers  is  also  presented  in  this  chapter.  Furthermore,  infor¬ 
mation  concerning  the  numerical  integration  and  time-history  analysis  performed 
by  the  computer  program  is  provided  at  the  end  of  the  chapter. 

Chapter  6  presents  the  most  important  results  of  all  the  analyses  undertaken 
during  this  research.  In  particular,  the  results  obtained  with  the  finite  element 
models  and  with  the  computer  program  developed  for  this  research  are  presented. 
A  comparison  of  the  results  is  also  included.  The  results  are  presented  in  the  form 
of  graphs  and  tables.  The  chapter  ends  with  a  discussion  of  the  results  obtained. 

Chapter  7  contains  the  conclusions  and  recommendations.  A  summary  of  the 
main  findings  and  achievements  is  discussed.  A  list  of  recommendations  is 
included,  to  provide  guidance  for  further  studies. 
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Building  Code  (International  Conference  of  Building  Officials,  ICBO  1997).  The 
majority  of  these  towers  are  very  lightly  reinforced  structures.  A  later  US  ACE 
study  (Dove  1998)  concentrated  efforts  on  the  evaluation  of  the  ductility  charac¬ 
teristics  of  these  towers.  The  final  goal  was  to  develop  nonlinear  analysis  tech¬ 
niques  to  estimate  the  ductility  and  to  determine  the  ability  of  these  towers  to 
resist  the  earthquake  demand.  An  additional  objective  was  to  identify  the  possible 
failure  mechanisms  in  these  towers  by  means  of  laboratory-scale  tests  (Dove 
2000). 

The  cost  of  retrofitting  an  existing  tower  has  been  estimated  to  be  between 
$5  million  and  $100  million  depending  on  the  work  required  (Dove  2000). 
Therefore,  an  accurate  estimation  of  the  seismic  performance  could  potentially 
lead  to  significant  savings.  Nevertheless,  the  retrofitting  is  an  expensive  process, 
and  it  causes  disruptions  to  the  operation  of  the  dam.  To  avoid  these  problems, 
intake  towers  in  seismically  active  regions  should  be  designed  to  withstand  earth¬ 
quakes  using  rational  analytical  methods  based  on  a  thorough  understanding  of 
the  dynamic  behavior  of  such  towers  (Daniell  and  Taylor  1994).  The  success  of  a 
given  tower  in  resisting  a  seismic  event  depends  on  the  magnitude  of  the  loads 
imposed  by  the  earthquake  and  the  structural  details  that  control  the  characteris¬ 
tics  of  the  nonlinear  dynamic  response.  It  has  been  argued  that  the  techniques 
currently  used  for  dynamic  analysis  of  intake  towers  may  not  properly  include 
these  factors  (Dove  1996). 


Seismic  Response  of  Intake-Outiet  Towers 

According  to  previous  studies,  the  response  of  an  intake  tower  under 
dynamic  loading  can  be  strongly  influenced  by  its  interaction  with  the  surround¬ 
ing  outside  and  inside  water  (Daniell  and  Taylor  1994).  It  is  now  accepted  that 
the  surrounding  water  modifies  the  dynamic  properties  of  the  tower  and  causes 
additional  dynamic  forces  (Chopra  and  Liaw  1975).  One  of  the  leading  authori¬ 
ties  on  the  subject.  Prof.  Anil  Chopra,  has  shown  that  the  water  acts  as  an 
incompressible  fluid  in  the  important  modes  of  vibration  of  the  tower.  The 
hydrodynamic  interaction  effects  can  then  be  represented  by  equivalent  added 
masses  along  the  height  of  the  tower.  These  effects  are  taken  into  account  exactly 
by  added-mass  functions  if  the  tower  is  uniform,  and  quite  accurately  if  the  cross- 
sectional  dimensions  of  the  tower  vary  along  its  height  (Goyal  and  Chopra 
1989b).  The  surrounding  water  tends  to  lower  the  vibration  frequencies  and  to 
modify  the  mode  shapes  (Chopra  and  Liaw  1975).  The  simple  approach  based  on 
added  masses  may  introduce  errors  at  higher  vibration  modes,  but  for  a  response 
dominated  by  the  fundamental  vibration  modes,  it  is  adequate  to  adopt  this  sim¬ 
plified  procedure.  Williams  (1990)  analyzed  the  seismic  response  of  intake-outlet 
towers  using  a  more  sophisticated  approach,  i.e..  Green’s  functions  to  obtain  the 
response  of  the  towers.  Williams  placed  special  emphasis  on  the  study  of  squatty 
towers  (towers  with  larger  diameter  to  water  depth  ratios).  In  his  investigation, 
Williams  demonstrated  that,  for  both  slender  and  squatty  towers,  the  interior 
surrounding  water  tends  to  lower  the  natural  frequencies  of  the  system  compared 
with  those  predicted  by  including  only  the  exterior  surrounding  water.  The  proper 
consideration  of  the  fluid-stmcture  interaction  effects  is  important  to  imderstand 
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the  behavior  of  the  intake  towers  during  the  event  of  an  earthquake.  Although 
important,  this  topic  is  not  the  primary  focus  of  this  research.  For  the  analyses  to 
be  undertaken,  die  hydrodynamic  interaction  will  be  considered  using  accepted 
modem  procedures. 

It  was  also  demonstrated  that  the  earthquake  response  of  intake  towers  can  be 
computed  with  enough  accuracy  by  considering  only  die  contribution  of  the  first 
natural  vibration  modes  (Goyal  and  Chopra  1989c).  Tliis  is  ^reeptable  when  their 
fundamental  vibration  periods  are  located  in  the  acceleration  or  velocity- 
controlled  regions  of  the  earthquake  response  spectrum.  In  another  sludy,  the 
same  authors  (Goyal  and  Chopra  1989b)  demonstrated  th^  the  influence  of  the 
shear  deformations  and  rotatory  inertia  should  be  considered  in  the  analysis  of 
squat  towers  when  only  the  first  two  vibration  modes  are  taken  into  account. 

Soil-stracture  interaction  effects  constitute  another  important  factor  that  may 
have  a  significant  influence  on  the  seismic  response  of  intice  towers.  Ihe  tower- 
foundation-soil  interaction  increases  the  fundamental  natural  period  and  the 
effective  damping  at  this  period.  Ihese  effects  are  the  consequence  of  the  flexi¬ 
bility  and  energy  dissipation  in  the  foimdation-soil  region  (Goyal  and  Chopra 
1989c).  ITie  dissipation  of  energy  through  radiation  (and  to  a  lesser  extent,  the 
material  damping)  h^  a  significant  effect  in  the  response.  The  tower-foundation 
interaction  effects  are  stronger  for  stiff  stractures  and  for  those  with  large  tower 
height  to  footing  radius  ratio.  The  influence  of  this  effect  depends  on  several 
factors  such  as  the  flexibility  of  tire  foundation  soil  and  the  slenderness  ratio,  and 
it  becomes  smaller  in  the  presence  of  water. 


Multiple-Support  Seismic  Excitation 

Free-standing  intake  towera  are  typically  connected  to  an  access  bridge,  and 
this  feature  may  influence  their  response  during  an  earthquake.  Ihis  response 
occurs  because  the  supports  of  the  tower  and  bridge  are  located  at  different  eleva¬ 
tions  and  thus  may  experience  different  groimd  accelerations.  Therefore,  to  <rarry 
out  a  more  rigorous  analysis  of  stractures,  it  is  necessary  to  corwider  the  effects 
of  the  spatially  varying  groimd  motions  at  the  supporte. 

The  effects  of  the  variation  of  the  support  motions  on  the  response  of 
extended  stractures  have  been  investigated  for  some  time.  The  principal  causes  of 
the  spatial  variation  in  seismic  ground  motions  are  (Der  Kiureghian  and  Neuen- 
hofer  1991)  the  difference  in  the  arrival  times  of  seismic  waves,  the  loss  of 
coherency  of  the  motions  due  to  reflections  and  refi'actions  of  the  waves  in  the 
medium  where  they  travel,  and  the  difference  in  the  local  soil  conditions  at  each 
support  location.  Observations  during  past  earthquakes,  for  example  the  Loma 
Prieta  earthquake  of  17  October  1989,  demonstrated  the  importance  of  the  effects 
of  multiple  excitations  in  the  supports  of  some  extended  stractures  such  as 
bridges.  The  feilure  of  several  extended  stractures  during  the  Loma  Prieta  earth¬ 
quake  highlighted  the  need  for  a  better  understanding  of  this  response  (Der 
Kiureghian  and  Neuenhofer  1991). 
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To  understand  the  response  of  structures  subjected  to  multiple-support  exci¬ 
tations,  it  is  necessary  to  introduce  the  equation  of  motion  for  these  types  of 
systems.  The  equation  of  motion  for  a  stmcture  subjected  to  multiple-support 
excitations  can  be  approximated  by  Equation  1  (Rutenberg  and  Heidebrecht 
1988): 


Mil  +  Cu  +  Ku  =  -MRiig  (t) 


(1) 


where 


M  =  mass  matrix  (lumped  or  consistent)  associated  with  active  (or 
dynamic)  degrees  of  freedom  (dofs) 

u  =  relative  displacement  vector 

C  =  damping  matrix 

K  =  stiffiiess  matrix 

R  =  quasi-static  influence  matrix  (given  by  Equation  2) 

Mg  =  seismic  acceleration  vector  of  support  points 

^  =  (2) 

In  Equation  2,  is  the  inverse  of  the  stiffiiess  matrix  related  to  the  motion  at 
the  interior  dofs,  and  is  the  stiffiiess  matrix  whose  rows  are  associated  with 
the  interior  dofs  and  its  columns  with  the  motions  of  the  supports. 

Several  methods  and  procedures  have  been  developed  to  accoimt  for  the 
multiple-support  excitation  effects  in  extended  stractures  such  as  bridges,  dams, 
and  pipeline  systems.  To  select  the  proper  approach  for  the  analysis  of  intake 
towers  subjected  to  multiple-support  excitation,  a  survey  of  the  available  proce¬ 
dures  was  performed.  The  most  relevant  methods  are  briefly  reviewed  in  the 
following  paragraphs. 

Various  approaches,  such  as  response  spectrum  and  time-history  methods, 
were  used  to  evaluate  the  response  of  multi-degree  of  freedom  (MDF)  systems  to 
multiple-support  excitations  (Y amamura  and  Tanaka  1990).  One  of  the  disadvan¬ 
tages  of  the  time-histoiy  method  is  that  the  results  from  the  analysis  are  specific 
to  the  set  of  selected  time-histories.  Also,  this  method  requires  extensive  compu¬ 
tational  effort  (Der  Kiureghian  and  Neuenhofer  1991).  However,  if  the  motions 
at  all  the  supports  of  the  structure  are  available  as  a  function  of  time,  this  method 
accoimts  for  the  effect  in  an  exact  way.  Moreover,  it  is  easy  to  understand,  and  its 
implementation  is  not  very  difficult. 

Another  method  used  for  the  analysis  of  systems  subjected  to  spatially  vary¬ 
ing  ground  motion  is  the  random  vibration  approach.  This  method  accounts  for 
the  effects  in  the  correlation  between  support  motions  by  means  of  the 
corresponding  power  spectral  density  and  correlation  functions.  This  approach  is 
based  on  the  statistical  characterization  of  a  set  of  ground  motions  between  the 
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different  supports  of  a  structural  system  (Der  Kiureghian  and  Neuenhofer  1991). 
However,  this  method  is  not  presently  used  by  practicing  engineers. 

Der  Kiureghian  and  Neuenhofer  (1991)  proposed  a  response  spectrum 
method  that  properly  includes  the  efifeets  of  the  variability  in  ground  motions  that 
arise  from  the  wave  passage  effeet,  the  ineoherence  effect,  and  the  local  soil 
conditions.  This  method  is  based  on  the  random  vibration  theory.  The  procedure 
requires  the  calculation  of  cross-mode  eonelation  coeffieients,  influence  and 
modal  participation  factors,  and  peak  ground  response  at  each  station  correspond¬ 
ing  to  each  modal  damping  and  frequency.  In  theory,  the  method  is  straight¬ 
forward  because  it  employs  a  simple  eombination  rale.  However,  it  requires 
extensive  eomputation  and  the  knowledge  of  many  factors  lhat,  in  the  majority  of 
the  cases,  are  unknown. 

Another  study  by  Harichandran  and  Wang  (1988),  confirmed  that  the 
response  of  long  structures  (such  as  bridges,  dams,  and  pipeline  systems)  to 
earthquake  ground  exeitation  is  sensitive  to  traveling  wave  effects  as  w^ell  as 
spatial  correlation  effects  of  ground  accelerations.  Die  same  authors  (Harichan¬ 
dran  and  Wang  1989)  examined  the  effects  of  the  spatial  variation  of  the  ground 
motion  on  an  indeterminate  structure  using  as  an  example  a  two-span  beam 
model.  This  study  also  showed  that  it  is  important  to  consider  the  spatial  vari¬ 
ation  of  earthquaJke  ground  motion  in  the  analysis  of  long,  statically  indetermin¬ 
ate  structures. 

Zerva  (1991)  studied  die  effects  of  the  spatial  variability  and  propagation  of 
seismic  groimd  motions  in  lifeline  systems.  He  indicated  that  the  variability  of 
ground  motions  between  supporte  may  be  attributed  to  the  loss  of  coherence  and 
phase  difference  of  die  seismic  motions.  Die  loss  of  coherence  may  result  from 
the  variability  in  the  rapture  process  and  discontinuities  in  the  travel  path  of  the 
seismic  waves  from  the  source  to  the  site  (Zerva  1991).  The  same  author  (Zerva 
1988)  had  previously  examined  the  vertical  response  of  a  two-span  bridge  due  to 
the  spatial  variation  of  ground  motion.  Die  analysis  was  made  using  the  powder 
spectral  density  of  the  acceleration  and  the  correlation  function.  This  study 
showed  that  the  assumption  of  equal  support  motion  in  die  seismic  evaluation  of 
bridges  may  underestimate  the  response.  Die  same  author  (Zerva  1992)  also 
presented  a  review^  of  the  techniques  used  for  the  modeling  of  the  spatial  varia¬ 
bility  of  the  seismic  ground  motions.  Zerva  described  in  this  paper  the  procedure 
for  the  stochastic  analysis  and  assigned  special  importance  to  the  following  tasks: 

a.  Die  definition  of  the  smoothed  cross  covariance  fimctions  of  die  motion 
between  two  stations. 

b.  The  evaluation  of  the  smoothed  cross  spectrum  by  die  Fourier  transform. 

c.  The  estimation  of  the  coherency  (a  complex  number)  and  the  coherence 
(a  real  number  between  0  and  1).  Die  coherency  is  a  measure  of  similar¬ 
ity  in  the  seismic  motions.  For  unconelated  motions,  its  vdue  is  equal  to 
zero. 
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d.  The  estimation  of  the  lagged  coherency,  which  describes  the  phase  varia¬ 
bility  in  the  acceleration  time  series. 

The  study  also  presented  several  empirical  coherency  models.  The  majority  of 
these  models  are  event-specific  and  time-consuming,  and  hence  their  use  is  lim¬ 
ited  in  earthquake  engineering.  For  these  reasons,  the  models  will  not  be  consid¬ 
ered  for  this  investigation. 

Rutenberg  and  Heidebrecht  (1988)  also  proposed  a  response  spectrum  analy¬ 
sis  for  multiple-support  excitations  in  bridges  using  a  traveling  wave  approach. 
Their  method  approximates  the  modal  displacement  time-history  plot  near  its 
peak  using  a  sine  curve.  It  was  observed  that  the  accuracy  of  the  results  obtained 
with  this  method  depends  on  the  transit  time  of  the  wave  and  the  natural  period. 

A  good  agreement  in  the  results  was  found  in  the  range  of  Tmax<  0.25T  and  T 
>01  sec,  where  imax  is  the  transit  time  and  T  is  the  lowest  natural  period  affecting 
the  displacement  response. 

Loh  and  Ku  (1995)  developed  another  method  to  calculate  the  structural 
response  for  multiple-support  seismic  excitations,  based  on  the  studies  of  Der 
Kiureghian  and  Neuenhofer  (1991).  Loh  and  Ku  (1995)  emphasized  the  impor¬ 
tance  of  the  spatial  variation  of  the  ground  motion  and  the  procedures  for  com¬ 
bining  the  responses.  This  method  simplifies  the  power  spectral  density  of  the 
ground  motion  and  system  transfer  function  in  order  to  evaluate  analytically  the 
frequency  integrals  that  define  the  cross-correlation  coefficients.  The  authors 
pointed  out  that  they  were  able  to  achieve  a  great  reduction  in  the  computational 
time  of  the  response.  Also,  they  mentioned  that  even  though  they  neglected  the 
cross-correlation  effects  between  quasi-static  and  dynamic  responses,  this  did  not 
introduced  significant  errors. 

Shaw  (1975)  studied  the  seismic  response  of  nuclear  plant  equipment 
accounting  for  multiple-support  excitations.  Pieces  of  equipment  such  as  pipe¬ 
lines  are  subjected  to  spatial  variation  in  the  acceleration  of  the  supports  during 
seismic  excitation.  Shaw  proposed  a  response  spectrum  analysis  based  in  the 
typical  equations  of  motion.  The  proposed  method  decomposes  the  multiple- 
support  excitation  problem  into  a  series  of  problems  subject  to  a  single  support 
excitation.  Afterwards,  the  individual  responses  are  combined  using  the  square 
root  of  the  sum  of  the  squares  (SRSS)  rule. 

Loh  (1991)  studied  the  spatial  variability  of  ground  motion  in  buried  pipe¬ 
lines.  In  this  study,  he  put  special  emphasis  on  the  differences  in  soil  stiffness  and 
the  variations  in  seismic  waves  along  the  segments  of  the  pipelines.  He  used  a 
two-dimensional  model  to  represent  the  pipeline  system  subjected  to  different 
seismic  waves  (Figure  3).  In  this  model  Kp  and  m  represent  the  stiffiiess  and  mass 
of  the  pipeline,  respectively,  and  Kgi  and  Kg2  are  the  stiffness  of  the  soil  at  each 
segment.  Loh  demonstrated  the  importance  of  considering  the  variability  in  soil 
stiffness  to  calculate  the  response  of  this  type  of  system.  He  indicated  that  this 
effect  is  more  significant  than  the  spatial  variation  in  the  seismic  waves. 

A  more  recent  study  by  Zembaty  and  Rutenberg  (2002)  focused  on  the  local 
site  amplification  effects  on  the  seismic  response  of  multi-support  structures.  The 
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Figure  3.  Two-dimensional  model  of  a  pipeline  system 
(adapted  from  Loh  1991) 

authors  recognized  the  importance  of  three  main  sources  of  spatial  variability:  the 
complexity  of  the  seismic  source,  the  velocity  of  propagation  of  waves,  and  the 
geological  and  geometrical  heterogeneities  of  the  ground.  The  effects  of  the  first 
two  sources  on  the  structural  response  increase  for  structures  with  larger  dimen¬ 
sions.  However,  the  third  source  can  be  important  at  short  distances.  The  local 
soil  conditions  amplify  the  structural  response  and  may  contribute  substantially 
to  the  overall  spatial  seismic  effects.  The  authors  cited  the  collapse  of  the 
Cypress  Street  Viaduct  during  the  1989  Loma  Prieta  earthquake  as  an  example  of 
spatial  variability  of  the  motion  at  the  supports.  It  is  now  believed  that  the  signif¬ 
icant  variation  in  local  soil  conditions  along  the  viaduct  generated  different  exci¬ 
tations  at  the  supports,  causing  its  collapse.  They  adopted  for  their  study  a 
complex  coherency  model  in  order  to  measure  the  similarity  between  motions. 
These  authors  proposed  a  simplified  response  spectrum  approach  that  considered 
the  local  variations  in  site  amplification  based  on  the  studies  by  Safak  (1995). 

This  model  properly  reflects  the  physics  of  soil  amplification  resulting  from  the 
reflections  of  the  shear  waves.  The  study  by  Zembaty  and  Rutenberg  (2002) 
demonstrated  the  importance  of  the  site  amplification  effects  in  the  response  of 
multi-supported  structures.  These  effects  significantly  modify  the  seismic 
response  of  such  structures.  Zembaty  and  Rutenberg’s  study  serves  as  additional 
motivation  and  justification  for  this  technical  report. 

After  examining  the  literature  on  the  subject,  it  was  decided  to  carry  out  a 
time-histoiy  analysis  of  the  tower  plus  bridge  system  in  which  the  different  exci¬ 
tations  at  their  two  supports  are  considered.  This  will  allow  us  to  assess  the 
importance  of  the  effect  for  the  particular  case  under  study.  Although  all  the 
previously  mentioned  investigations  concluded  that  the  influence  of  the  multiple- 
support  excitations  on  the  response  is  important,  they  dealt  primarily  with  spa¬ 
tially  extended  structures.  The  tower  plus  bridge  system  is  not  quite  a  spatially 
extended  system;  the  difference  in  the  support  motions  will  arise  from  the  fact 
that  the  supports  are  at  different  elevations  along  the  dam.  As  shown  in  Figure  2, 
the  base  of  the  intake  tower  is  near  the  toe  of  the  dam,  and  the  abutment  of  the 
bridge  is  close  to  the  crest. 
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3  General  Methodology 


An  overview  of  the  general  methodology  and  a  brief  description  of  the  finite 
element  models  used  during  this  research  are  presented  in  this  chapter.  Several 
analytical  tools  were  used  to  accomplish  the  objectives.  The  study  began  with  a 
review  of  the  methods  currently  available  for  the  analysis  of  structures  with 
multiple-support  excitations.  This  study  was  done  to  find  the  most  practical 
method  for  the  analysis  of  the  system  formed  by  the  intake  tower  and  the  bridge. 
A  review  of  the  studies  related  to  intake-outlet  towers  and  access  structures  was 
also  performed.  This  study  was  used  to  obtain  general  knowledge  about  the  fac¬ 
tors  that  affect  the  seismic  response  of  this  structural  system.  These  factors 
included  the  hydrodynamic  and  soil-structure  interaction  effects.  This  review  was 
also  intended  to  identify  methods  for  the  analysis  and  computer  modeling  of 
these  towers.  A  typical  existing  intake-outlet  tower  was  selected  for  the  case 
study. 

The  next  step  of  the  research  was  the  development  of  several  finite  element 
models.  The  finite  element  models  were  created  using  the  SAP2000  (CSI 2000) 
and  QUAD4M  (Idriss  et  al.  1973)  computer  programs.  These  models  included 

a.  A  frame  model  of  the  intake  tower  using  SAP2000. 

b.  A  finite-element  model  of  the  earthen  dam  using  QUAD4M. 

c.  A  finite-element  model  of  the  earthen  dam  using  SAP2000. 

d.  A  finite-element  model  of  the  combined  dam-tower  system  using 
SAP2000. 

e.  A  frame  model  of  the  intake  tower  with  added  springs  and  dashpots 
representing  soil-structure  effects  using  SAP2000. 

The  finite  element  models  were  used  for  the  seismic  response  analysis  of  the 
intake  tower  and  bridge  system.  An  acceleration  time-history  was  applied  to  the 
base  of  each  finite  element  model.  The  same  boundary  conditions  were  applied  to 
each  model  to  compare  the  results.  These  models  were  developed  to  obtain  the 
acceleration  time-histories  at  the  base  of  the  tower  and  at  the  support  of  the 
inspection  or  access  bridge.  Detailed  information  about  the  conditions  used  in 
each  model  will  be  discussed  in  the  next  chapter.  Only  a  brief  description  of  the 
models  is  presented  in  this  chapter  as  part  of  the  general  methodology  used. 
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Another  tool  developed  during  this  investigation  was  a  computer  program  for 
the  analysis  of  intake-outlet  towers  created  using  the  MATLAB  (MathWorks 
2000)  computational  environment.  The  program  can  account  for  the  effects  of 
multiple-support  excitations  by  means  of  a  time-history  analysis.  It  reads  the 
earthquake  accelero^ams  at  two  different  locations  in  the  horizontal  and  vertical 
directions.  The  formulation  for  the  response  calculation  is  based  on  the  relative 
displacements  formulation  of  the  equations  of  motion.  The  program  also  includes 
the  water-tower  interaction  effects  using  the  formulation  developed  by  Goyal  and 
Chopra  (198%),  which  introduces  the  concept  of  added  hydrodynamic  masses 
along  tile  height  of  the  tower.  The  program  can  include  the  access  structure  using 
truss  or  frame  elements.  In  the  last  phase  of  the  research,  the  soil-stiucture  inter¬ 
action  effects  were  included  in  this  computer  program.  The  problem  of  soil- 
structure  interaction  was  introduced  by  adding  the  springs  and  d^hpots  at  the 
base  of  the  tower  and  at  the  support  of  the  bridge. 

The  results  obtained  with  the  MATLAB  computer  program  and  with  the 
more  comprehensive  finite  element  models  were  examined  and  compared.  The 
absolute  acceleration  and  relative  and  absolute  displacement  time-histories  were 
compared.  Conclusions  and  recommendations  for  the  seismic  analysis  of  the 
towers,  including  the  importance  of  considering  the  bridge  and  the  multiple- 
support  excitations,  were  drawn  based  on  the  results. 


Selection  of  the  Case  Study 

A  specific  case  study  was  needed  as  framework  for  this  research  effort,  and 
the  earthen  dam  selected  was  Arkabutla  Dam.  This  example  was  selected  based 
on  recommendations  from  personnel  at  the  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC)  during  a  summer  internship  that  the  primary  author 
completed  at  the  ERDC  Geotechnical  and  Structures  Laboratory  (GSL)  in  Vicks¬ 
burg,  MS. 

Arkabutla  Dam  is  located  in  northern  Mississippi,  30  miles'  south  of 
Memphis,  TN,  west  of  Interstate  55  (Figure  4).  The  dam  is  of  a  rolled  earth  fill, 
approximately  1 1,500  ft  long  with  an  average  height  of  66  ft.  It  was  constructed 
during  the  period  1940  to  1943.  The  dam  is  founded  in  a  valley  filled  with  allu¬ 
vium.  The  alluvium  consists  of  a  top  stratum  of  fine-grained  sediments  and  a 
substratum  of  predominantly  coarse-grained  sediments  (Chang  1990).  The  fine¬ 
grained  sediments  consist  of  clays,  silts,  and  some  fine  sands. 

The  soil  properties  vary  along  the  length  of  the  dam.  The  soil  strata  were 
divided  into  17  layers  to  obtain  a  better  representation  of  the  dam  in  the  finite 
element  model.  These  properties  were  obtained  from  USACE  data  available  in 
the  GSL.  Table  1  presents  the  values  for  the  shear  modulus  (G),  shear  wave 
velocity  (Vs),  unit  weights  (y),  and  the  number  of  shear  and  damping  curves  used 
in  the  analysis  for  the  different  types  of  soils.  The  modulus  and  damping  degra¬ 
dation  curves  are  presented  in  Figures  5  and  6,  respectively.  These  plots  show, 

'  A  table  of  factore  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  viii. 
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Figure  4.  Location  of  Arkabutia  Dam 


respectively,  how  the  modulus  decreases  and  the  damping  ratio  increases,  due  to 
the  nonlinear  behavior  of  the  soil.  Theses  plots  were  obtained  from  the  values 
shown  in  Table  2.  This  table  shows  the  values  of  ratio  G/Gmax  and  damping  ratio 
as  a  function  of  the  shear  strain  y  used  as  input  to  the  model  developed  in 
QUAD4M  for  all  the  types  of  soil  that  comprise  the  dam. 

The  finite  element  model  of  the  embankment  was  developed  using  these  soil 
properties  in  two  computer  programs,  SAP2000  and  QUAD4M.  Figure  7  shows 
the  soil  strata  in  the  finite  element  model  of  the  Arkabutia  Dam  identified  accord¬ 
ing  to  the  material  number  listed  in  Table  1 . 

Initially,  the  real  intake-outlet  tower  of  the  Arkabutia  Dam  was  chosen  for 
the  analysis.  However,  because  the  geometry  of  this  tower  is  quite  complex,  a 
typical  intake-outlet  tower  representative  of  the  existing  tower  inventory  in  the 
United  States  was  finally  used.  In  addition,  using  this  typical  tower  allows  a 
general  understanding  of  the  dynamic  behavior  of  the  whole  population  of 
USAGE  towers.  The  dimensions  were  chosen  based  on  data  available  from  the 
Corps’  inventory  (Figure  8). 

Dove  (1996)  conducted  a  comprehensive  study  aimed  at  identifying  the  most 
important  structural  characteristics  of  intake  towers  in  the  USAGE  inventory. 

This  study  focused  on  162  intake  towers.  A  summary  of  the  results  is  presented 
in  Table  3,  highlighting  the  mean  value  and  standard  deviation  of  some  of  the 
main  geometric  properties  of  these  towers  (Dove  and  Matheu  2000a,b).  Based  on 
these  dimensions,  typical  tower  and  access  structures  were  defined  to  use  in  con¬ 
junction  with  the  analyses  of  Arkabutia  Dam.  Table  4  summarizes  the  final 
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Table  2  (Continued) 


Curve  1(3  Modulus  Reduction  SAND  Average  (Seed  and  Idriss  1970) 


0.0001 

0.0003 

0.001 

0.003 

0.01 

0.03 

0.01 

0.3 

1 

1 

0.98 

0.95 

0.89 

0.73 

0.52 

0.29 

0.14 

0.06 

Curve  #3  Damping _ SAND  Average  (Seed  and  Idriss  1970) 


Strain,  y  (%) 
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Curve  #4 

Modulus  Reduction  SAND  Average  (Seed  and  Idriss  1970) 
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Curve  #4 

Damping 

SAND  Average  (Seed  and  Idriss  1970) 
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Curve  #5 

Modulus  Reduction  CLAY  (Pl=20  Sun  et  al.  1988) 
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Curve  #5 

Damping 

CLAY  Lower  Bound  (Seed  and  Idriss  1970) 
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Modulus  Reduction  CLAY  (PI=20-40  Sun  et  al.  1988) 
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Curve  #6 

Damping 

CLAY  Average  (Seed  and  Idriss  1970) 
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Curve  #7 

Modulus  Reduction  CLAY  (FI=20-40  Sun  et  al.  1988) 
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Damping 
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Table  2  (Continued) 


Curve  #8  Modulus  Reduction  CLAY  (Pl=20-40  Sun  et  al.  1970) 


strain,  y  (%) 
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Curve  #8 _ Damping _ SAND  Average  (Seed  and  Idriss  1970) 
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Curve  #9 

Modulus  Reduction  CLAY  (PI=:20-40  Sun  et  al.  1988) 
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Damping 
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Curve  #10 

Modulus  Reduction  CLAY  (PN20-40  Sun  et  al.  1988) 
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Damping 

CLAY  Average  (Seed  and  Idriss  1970) 
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Curve  #1 1 

Modulus  Reduction  CLAY  (PI=20-40  Sun  et  al.  1988) 
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Curve  #11 

Damping 

CLAY  Average  (Seed  and  Idriss  1970) 
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Curve  #12 

Modulus  Reduction  CLAY  (Pl=20-40  Sun  et  al.  1988) 
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Figure  8.  Dimensions  of  Arkabutia  Dam 


Table  3 

Parameters  of  the  Intake  Towers  of  the  USAGE  Inventory 

1  Parameter 

Mean  Value 

Std  Deviation 

Total  structural  height,  ft 

163.31 

19,31 

Section  width  (parallel  to  flow),  ft 

32.86 

3.91 

Section  width  (perpendicular  to  flow),  ft 

35.69 

3.88 

Thickness  of  the  walls  (parallel  to  flow),  ft 

3.21 

0.64  1 

Thickness  of  the  walls  (perpendicular  to  flow),  ft 

3.28 

0.63  1 

Vertical  steel  reinforcement  (walls  parallel  to  flow),  % 

0.273 

0.178 

Vertical  steel  reinforcement  (walls  perpendicular  to  flow),% 

0.275 

0.166 

Table  4 

Final  Parameters  of  the  Typical  Intake-Outlet  Tower 

Parameter _  Value 

Total  staictural  height,  ft  75.00 

Section  width  (parallel  to  flow),  ft  30.00 

Section  wdth  (perpendicular  to  flow),  ft  30.00 

Thickness  of  the  walls,  ft  3.00 

Height  of  the  outside  water,  ft  57.00 

Height  of  the  inside  water,  ft  57.00 

Height  of  the  crest  dam,  ft  66.00 


parameters  and  dimensions  selected  for  the  case  study  used  in  this  research.  The 
final  access  bridge  selected  was  longer  than  the  actual  bridge  of  the  dam:  it  has 
an  approximate  length  of 202.0  ft.  The  actual  access  bridge  is  a  truss-type 
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structure;  an  equivalent  rectangular  section  was  chosen  for  the  bridge  in  order  to 
simplify  the  analysis  and  to  generate  the  computer  models  in  two  dimensions. 
The  intake-outlet  tower  selected  has  a  prismatic  hollow  square  section  along  its 
height.  The  materials  of  the  tower  and  service  bridge  are  concrete  and  steel, 
respectively.  Figure  9  shows  the  cross  section  and  front  view  of  the  tower 
selected. 


Typical  Hake  Tower 


Figure  9.  Cross  section  and  front  view  of  the  typical  intake-outlet  tower  (units  in  feet) 


As  mentioned  previously,  this  tower  was  used  in  all  the  finite  element  and 
frame  models  developed  for  this  investigation.  It  is  expected  that  this  case  study 
will  allow  a  better  understanding  of  the  dynamic  behavior  of  intake-outlet  towers 
with  similar  characteristics.  This  analysis  also  intends  to  identify  the  factors  that 
have  more  influence  on  the  seismic  response  of  this  type  of  structural  system. 


Summary 

This  chapter  presented  a  brief  description  of  the  methodology  used  in  this 
research.  Basically,  this  methodology  can  be  divided  into  four  phases;  literature 
review,  definition  and  development  of  the  finite  element  models,  generation  of 
the  computer  program  using  the  MATLAB  computational  environment,  and  anal¬ 
ysis  of  the  results.  In  the  next  chapter,  a  detailed  description  of  the  conditions 
used  in  each  finite  element  model  will  be  presented. 
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The  selection  of  the  particular  case  study  for  this  investigation  was  also 
presented  in  this  chapter.  The  intake-outlet  tower  selected  was  based  on  typical 
dimensions  found  in  an  inventory  of  the  towers  under  USAGE  jurisdiction.  The 
dam  selected  for  this  study  was  Arkabutla  Dam  in  the  state  of  Mississippi.  This 
dam  was  chosen  because  of  the  interest  of  the  Corps  in  its  study.  Special  consid¬ 
erations  were  taken  to  adjust  the  typical  tower  and  access  structure  to  the  dimen¬ 
sions  of  Arkabutla  Dam. 
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4  Finite  Eiement  Models 


This  chapter  presents  a  description  of  the  finite  element  models  developed 
for  this  investigation.  The  parameters,  boundary  conditions,  type  of  analyses,  and 
input  data  used  in  each  model  are  also  described.  The  finite  element  models  as 
described  in  the  previous  chapter  included  (a)  a  frame  model  of  the  stand-alone 
intake-outlet  tower,  (b)  a  model  of  the  Arkabutla  Dam  using  the  program 
QUAD4M,  (c)  a  model  of  the  Arkabutla  Dam  using  the  program  SAP2000,  (d)  a 
model  of  the  Arkabutla  Dam  plus  the  intake-outlet  tower  using  SAP2000,  and 
(e)  a  spring-dashpot  model  of  the  tower  in  SAP2000  and  MATLAB.  Each  model 
was  analyzed  to  identify  the  importance  of  the  multiple-support  excitation  effects 
and  other  factors  in  the  response  of  these  towers.  Moreover,  descriptions  of  the 
computer  programs  SAP2000  and  QUAD4M  are  presented  for  informative 
purposes. 


Brief  Description  of  the  Computer 
Program  SAP2000 

Several  commercial  software  packages  have  been  developed  for  the  analysis 
and  design  of  all  types  of  civil  engineering  structures.  The  Structural  Analysis 
Program  (SAP)  is  one  of  them.  This  program  was  originally  developed  by 
Prof  Edward  Wilson  at  the  University  of  California  at  Berkeley  in  the  1970s.  It 
was  later  updated  and  commercialized  by  Computers  and  Structures,  Inc.,  of 
Berkeley,  CA.  The  most  recent  version,  SAP2000,  is  among  the  most  powerful 
and  versatile  programs  available  on  the  market.  Some  of  the  program  features 
include  static  and  dynamic  analysis,  linear  and  nonlinear  analysis,  dynamic  seis¬ 
mic  analysis,  static  pushover  and  P-Delta  analysis,  bridge  analysis,  nonlinear 
links,  spring  elements,  multiple  coordinate  systems,  many  types  of  constraints, 
and  several  other  features  (CSI 2000).  The  majority  of  the  models  can  be  created 
using  the  SAP2000  graphical  user  interface.  This  graphical  interface  makes  the 
program  very  user  friendly.  For  all  these  reasons,  this  computer  program  was 
selected  as  the  main  tool  for  the  analyses  performed  during  this  research. 
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Brief  Description  of  the  Computer 
Program  QUAD4M 

The  program  QUAD4M  was  used  during  this  research  to  develop  one  of  the 
finite  element  models  of  Arkabutla  Dam.  QUAD4M  is  a  modification  of  the 
QUAD4  program,  which  was  developed  in  1973  by  researchers  at  the  Depart¬ 
ment  of  Civil  Engineering,  University  of  California  at  Berkeley  (Idriss  et  al. 
1973).  The  computer  programs  QUAD4  and  QUAD4Mwere  developed  to  evalu¬ 
ate  the  seismic  response  of  soil  structures  using  finite  element  procedures.  The 
modifications  included  in  the  QUAD4M  were  the  implementation  of  trammitting 
boundaries  and  the  improvement  of  the  time-stepping  algorithm  and  seismic 
coefficient  calculations  (Hudson  et  al.  1994).  The  formulation  of  the  program 
includes  a  variable  damping  finite  element  procedure.  This  feature  permits  the 
accurate  analysis  of  soil  structures  with  irregular  boundary  conditions  and  vari¬ 
ations  in  material  properties.  The  irregular  boundary  conditions,  when  present, 
cause  the  strains  developed  in  each  element  to  vary  considerably  compared  with 
one-dimensional  models  such  as  those  used  in  the  SHAKE  program.  Because  the 
damping  in  soils  is  strain  dependent,  the  damping  value  should  be  based  on  the 
strain  developed  in  each  element.  Thus,  the  use  of  a  constant  damping  value  can 
lead  to  inaccurate  results  for  some  of  those  conditions.  In  addition,  this  procedure 
allows  for  the  incorporation  of  both  stiffiiess  and  damping  values  that  are  strain 
dependent,  for  each  element  (Idriss  et  al.  1973).  QUAD4M  considers  the  nonlin¬ 
ear  behavior  of  the  soils  in  an  approximate  way  using  the  Equivalent  Linear 
Method.  The  output  from  the  program  includes  acceleration,  displacement,  and 
shear  strain  time-histories.  Because  the  program  is  very  powerful  for  analyzing 
soil  structures,  it  was  used  to  develop  the  first  finite  element  model  of  Arkabutla 
Dam. 


Finite  Eiement  Models 

Five  analysis  cases  were  considered  in  this  research  to  evaluate  the  effects  of 
the  multiple-support  excitations  in  these  towers.  Other  factors  were  also  consid¬ 
ered  to  obtain  a  complete  panorama  of  the  dynamic  behavior  of  this  structural 
system.  These  factors  included  the  hydrodynamic  and  soil-structure  interaction 
effects. 

First,  a  simplified  model  of  the  intake-outlet  tower  was  developed.  This 
model  is  very  simple  and  consists  of  only  the  tower  itself  modeled  using  fi:ame 
elements.  In  the  current  engineering  practice,  similar  column  models  are  used  to 
study  intake  towers  when  their  geometry  is  relatively  simple.  This  model  was 
done  only  for  preliminary  analysis;  it  provides  an  overview  of  file  response  of  the 
tower  without  the  effects  of  the  access  bridge.  Moreover,  the  model  was  also 
used  to  assess  the  importance  of  the  inclusion  of  the  access  bridge  in  the  analysis 
of  the  towers  and  the  effects  of  the  spatial  variability  of  the  ground  motions  at  the 
supports.  The  model  of  the  tower  was  created  using  SAP2000. 

The  simplified  model  was  developed  according  to  the  geometry  described  in 
Chapter  3,  corresponding  to  a  prismatic  hollow  section  along  the  height  of  the 
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tower  with  top  and  base  slabs.  The  two-dimensional  model  (ZX  plane)  was  dis¬ 
cretized  using  a  total  of  12  frame  elements  and  13  nodes.  The  boundary  condi¬ 
tions  imposed  at  the  end  of  the  tower  were  fixed  in  three  directions  (translations 
in  X  and  Z,  and  rotation  around  Y)  whereas  the  other  points  of  the  tower  were 
free  to  move.  Figure  10  shows  the  resulting  simplified  model  obtained  from 
SAP2000. 


Figure  10.  Simplified  model  of  the  intake-outlet  tower  using 
SAP2000 


Other  considerations  were  taken  into  account  for  the  development  of  this 
model.  For  example,  shear  deformations  and  rotational  inertia  were  neglected  in 
the  analysis.  These  effects  were  neglected  in  order  to  compare  the  results  with  the 
MATLAB-based  program  developed  for  this  investigation.  The  formulation  of 
the  MATLAB  program  is  based  in  Bernoulli  beam  elements:  these  are  the  sim¬ 
plest  beam  elements  and  do  not  take  into  account  the  effects  of  the  shear  defor¬ 
mation  and  rotational  inertia.  The  beam  elements  in  SAP2000  are  based  in  the 
Timoshenko  beam  theory,  and  thus  the  frame  and  beam  elements  include  these 
two  special  effects. 

Fluid-structure  interaction  effects  were  included  in  the  model  using  the  con¬ 
cept  of  hydrodynamic  added  masses.  It  was  assumed  that  the  water  acts  as  an 
incompressible  fluid.  The  simplified  procedure  developed  by  Goyal  and  Chopra 
(1989a)  was  used  to  obtain  these  added  masses.  A  brief  description  of  this  proce¬ 
dure  is  presented  in  Chapter  5.  The  input  used  in  the  program  was  an  acceleration 
time-history  in  the  horizontal  direction.  Vertical  acceleration  time-histories  were 
not  considered  in  the  analyses  done  with  this  simple  model.  The  effects  of  the 
vertical  acceleration  in  intake-outlet  towers  are  neglected  in  the  majority  of  the 


Chapter  4  Finite  Element  Models 


27 


cases,  and  this  is  justified  because  these  structures  are  extremely  rigid  in  axial 
direction. 

To  evaluate  the  response  of  the  towei^  subjected  to  seismic  excitation,  a 
time-history  analysis  was  done.  The  relative  and  absolute  displacement  time- 
histories  were  obtained  from  the  analysis.  First,  a  modal  analysis  of  the 
simplified  model  was  performed.  From  the  modal  analysis,  the  mode  shapes, 
modal  periods,  and  frequencies  of  the  tower  were  obtained. 

The  principal  objective  of  this  research  was  to  examine  and  establish  the 
importance  of  the  effects  of  the  multiple-support  excitations  in  the  intake-outlet 
towers.  For  this  reason,  a  finite  element  model  of  Arkabutla  Dam  was  created.  To 
obtain  reliable  results,  it  is  important  to  choose  a  tool  that  is  developed  specifi¬ 
cally  for  soil  structures  such  as  earthen  dams  and  that  can  provide  accurate 
results.  Therefore,  the  program  QUAD4M  was  used  to  develop  the  first  model  of 
Arkabutla  Dam.  QUAD4M  analyzes  soil  structures  in  a  more  refined  fashion  than 
other  general-purpose  finite  element  programs,  and  it  also  includes  a  variable 
damping  finite  element  procedure.  This  model  was  used  to  calibrate  the  third 
model  developed  with  SAP2000.  Figure  1 1  shows  the  mesh  of  the  finite  element 
model  of  the  dam  developed  in  QUAD4M. 


Figure  11.  Finite  element  model  of  Arkabutla  Dam  using  QUAD4M 


The  finite  element  model  of  the  dam  consists  of  1,404  elements  and  1,479 
nodes.  Plane  strain  elements  were  used  in  the  model.  The  boundary  conditions 
imposed  at  the  base  of  the  embankment  were  fixed  in  the  two  directions  (transla¬ 
tion  in  the  horizontal  X  axis  and  in  the  vertical  Z  axis).  At  the  side  nodes  of  the 
dam  foundation,  only  the  vertical  translation  (Z  direction)  was  restricted  for  the 
case  of  horizontal  seismic  motion.  This  condition  was  achieved  with  rollers  in  the 
horizontal  direction.  The  model  was  divided  into  17  zones  representing  the  dif¬ 
ferent  soil  layers  (see  Chapter  3).  The  input  used  for  the  program  was  a  site- 
specific  acceleration  time-history. 

This  model  was  developed  to  obtain  the  absolute  accelerations  at  the  top  of 
the  dam  (point  of  contact  of  the  access  bridge)  and  at  die  base  of  the  tower.  The 
dynamic  analysis  starts  with  an  eigenvalues  analysis  (modal  analysis)  of  the  dam 
for  the  determination  of  the  modal  periods,  frequencies,  and  modes.  Then,  a 
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time-history  analysis  is  performed  to  calculate  the  acceleration  time-histories  at 
the  required  locations. 

The  third  finite  element  model  created  also  corresponds  to  Arkabutla  Dam, 
but  now  uses  SAP2000.  The  geometry  of  the  model  is  identical  to  the  model 
developed  with  QUAD4M.  The  creation  of  another  identical  finite  element 
model  was  needed  because  QUAD4M  cannot  consider  frame  structures  such  as 
the  intake  towers.  As  was  mentioned  before,  the  program  was  developed  to 
evaluate  and  analyze  soil  structures.  SAP2000,  on  the  other  hand,  can  combine 
plane  strain  elements  (such  as  the  dam)  with  frame  elements  (such  as  the  intake- 
outlet  tower).  The  first  analyses  of  the  dam  using  QUAD4M  were  done  because 
the  program  is  able  to  provide  a  better  approximation  of  the  seismic  response  of 
earth  structures.  This  situation  occurs,  in  part,  because  of  the  model’s  considera¬ 
tion  of  the  nonlinear  behavior  of  the  soil  and  the  use  of  absorbing  boundaries. 
Therefore,  the  SAP2000  model  was  calibrated  to  obtain  acceleration  time- 
histories  similar  to  the  ones  obtained  with  QUAD4M.  The  calibration  was 
accomplished  by  adjusting  the  value  of  the  modal  damping  parameter.  A  damp¬ 
ing  ratio  of  24.5  percent  was  obtained  from  the  calibration.  This  value  was  used 
for  all  the  modes  included  in  the  time-history  analysis  performed  with  SAP2000 
because  this  program  does  not  permit  assigning  different  damping  values  to  the 
elements.  Figure  12  illustrates  the  finite  element  model  generated  with  SAP2000, 


Figure  12.  Finite  element  model  of  Arkabutla  Dam  using  SAP2000 


The  model  was  done  in  Z-X  plane  (2-D).  The  same  boundary  conditions  as 
previously  described  for  the  QUAD4M  model  were  applied  to  the  base  of  the 
embankment.  The  nodal  points  at  the  base  of  the  foundation  were  fixed  both  in 
the  vertical  and  horizontal  directions.  Similar  to  the  other  model,  horizontal 
rollers  were  attached  to  the  nodes  at  the  two  sides.  Here  again,  only  horizontal 
acceleration  time-histories  were  used  as  input. 

The  two  models  of  the  dam  just  described  permit  the  study  of  the  dynamic 
behavior  of  the  dam  alone.  To  evaluate  and  investigate  the  response  of  the  dam 
with  the  intake-outlet  tower,  it  is  necessary  to  use  a  finite  element  model  that 
contains  both  structures.  This  combined  model  can  be  used  to  analyze  the  effects 
of  having  different  excitations  at  the  supports,  which  is  the  main  objective  of  the 
research.  In  addition,  the  tower  plus  dam  model  is  useful  to  consider  the 
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soil-structure  interaction.  The  mesh,  boundary  conditions,  materials,  md  input 
accelerations  used  in  this  model  were  the  same  as  those  considered  previously  for 
the  dam  alone.  The  only  modification  was  the  inclusion  of  the  intake  tower.  Fig¬ 
ure  13  illustrates  the  finite  element  model  of  the  Arkabutla  Dam  plus  the  intake 
tower. 


Figure  13.  Finite  element  model  of  Ari<abutla  Dam  plus  its  intake  using  SAP2000 


The  properties  of  the  tower  were  described  in  Chapter  3.  The  tower  was 
divided  into  12  frame  elements  and  13  nodes.  The  access  bridge  was  modeled  as 
a  long  truss  element,  whose  properties  were  presented  in  Chapter  3.  End  releases 
for  the  moments  were  imposed  to  join  the  frame  element  in  the  access  bridge  to 
the  truss  element.  The  access  bridge  is  connected  at  the  crest  of  the  dam;  its  ele¬ 
vation  from  the  base  of  die  tower  is  66  ft.  A  massless  and  very  rigid  frame  ele¬ 
ment  was  embedded  into  the  dam  from  the  base  of  the  tower  in  order  to  model 
more  accurately  the  fixed  base  condition  between  the  dam  and  tower. 

The  output  obtained  from  this  model  essentially  consists  of  the  absolute 
acceleration,  relative  and  absolute  displacement  time-histories,  and  modal  peri¬ 
ods  and  frequencies  of  the  combined  system. 

A  model  of  the  intake  tower  consisting  of  a  frame  structure  with  springs  and 
dashpots  at  the  base  was  developed  in  the  last  phase  of  this  research.  This  model 
was  used  to  take  into  account  the  soil-structure  interaction  effects  in  a  simpler 
way  than  in  the  tower-dam  model.  This  phenomenon  is  caused  by  the  contact 
stresses  between  soil  and  foundation,  which  are  responsible  for  the  deformations 
in  the  soil  (Newmark  and  Rosenblueth  1971).  Generally,  the  soil-structure  inter¬ 
action  has  little  effect  on  the  dynamic  response  of  many  structures  and  founda¬ 
tion  systems.  In  other  cases,  particularly  when  massive  structures  are  founded  in 
soft  soils,  this  interaction  effect  can  be  significant  (Kramer  1996).  Initially,  it  was 
not  planned  to  consider  this  effect  in  the  study.  However,  when  the  results  of  the 
fixed-base  model  created  in  MATLAB  were  compared  with  the  SAP2000  finite 
element  model,  there  were  discrepancies  that  could  be  attributed  only  to  the  soil- 
structure  interaction.  A  detailed  and  comprehensive  evaluation  of  the  results 
obtained  with  this  model  is  discussed  in  Chapter  7.  Figure  14  shows  the  frame 
with  springs-dashpots  model  of  the  intake  tower  that  was  created  to  approxi¬ 
mately  account  for  the  soil-structure  interaction. 
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Linear  springs  (which  represent  the  flexibility  of  the  soil)  were  attached  to 
the  base  of  the  tower  and  at  the  end  of  the  access  bridge  in  three  degrees  of  free¬ 
dom  (horizontal  and  vertical  translation  and  rocking).  Viscous  dashpots  were 
attached  at  the  base  of  the  tower  in  the  horizontal  and  vertical  directions.  Another 
dashpot  in  the  direction  of  the  rocking  motion  was  also  inserted  at  the  base.  No 
dashpots  were  placed  at  the  base  of  the  bridge.  This  was  done  because  there  is  no 
information  in  the  literature  to  define  the  damping  constants  of  the  dashpots  that 
represent  the  radiation  damping  at  the  crest  of  a  dam.  All  the  dashpot  coefficients 
available  in  papers  and  textbooks  were  obtained  for  a  uniform  or  layered  half¬ 
space. 

The  stiffness  and  damping  constants  for  the  elements  at  the  base  of  the  tower 
were  calculated  from  the  solutions  available  in  the  literature.  The  first  systematic 
solution  of  the  problem  of  defining  these  values  was  obtained  by  Lysmer  (1965). 
He  developed  the  stiffiiess  and  damping  constants  for  the  vertical  motion  of  a 
rigid  circular  plate  on  a  half-space  taking  into  account  the  proper  distribution  of 
contact  stresses.  Later  on,  other  researchers  proposed  similar  coefficients  for 
other  types  of  motion,  and  for  embedded  foundations,  foundations  with  other 
shapes,  and  foundations  on  layered  soil  deposits.  The  expressions  for  the  stiffiiess 
and  damping  coefficients  used  in  this  report  were  taken  from  the  technical  report 
for  SAP2000  (CSI 2000).  These  solutions  were  adapted  from  those  published  in 
the  book  Fundamentals  of  Earthquake  Engineering  (Newmark  and  Rosenblueth 
1971).  Table  10  of  the  book  displays  the  stiffness  and  damping  coefficients  for 
the  spring-dashpot  model  used  in  this  research.  These  constants  are  valid  for  a 
rigid  circular  plate  lying  on  the  surface  of  a  half-space  and  are  independent  of  the 
frequency  of  vibration. 
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In  Table  5,  G  is  the  small  strain  shear  modulus  of  the  soil,  r  represents  the 
plate  radius,  and  vand  p  are  the  Poisson’s  ratio  and  mass  density  of  the  soil, 
respectively.  When  a  noncircular  foundation  is  considered,  an  equivalent  radius 
must  be  defined  in  order  to  use  these  equations.  In  the  present  study,  the  equiva¬ 
lent  radius  was  obtained  by  equating  the  area  of  a  circular  plate  to  the  square 
plate  and  solving  for  r.  These  constants  were  introduced  to  the  spring-dashpot 
model  developed  in  SAP2000  and  in  MATLAB. 


Table  5 

Properties  of  Rigid  Circuiar  Plate  on  Surface  of  Half-Space 

Direction 

Stiffness 

Damping 

Vertical 

1-v 

C  =  1.79^iQcr" 

Horizontal 

(l-v*) 

^-18.2C?r(2_^y 

C  =  L0B^Kpr^ 

Rocking 

K^lJGr^ 

C= 0.47^ Kpr^ 

As  it  was  done  previously  for  the  other  cases,  a  time-history  and  modal  anal¬ 
ysis  were  performed  using  this  model.  Absolute  acceleration  and  absolute  and 
relative  displacement  time-histories  were  obtained.  Special  attention  was  paid  to 
the  time-histories  at  the  top  of  the  tower  to  compare  the  results  with  those  of  the 
full  tower  plus  dam  model. 


Summary 

This  chapter  presented  a  description  of  the  finite  element  models  developed 
for  this  investigation.  Five  models  v^ere  presented^  with  a  description  of  the 
boundary  conditions,  properties,  and  other  features  of  each  model.  For  all  the 
models,  a  time-histoiy  and  modal  (eigenvalue)  analysis  were  done.  These  analy¬ 
ses  allow  one  to  understand  the  seismic  response  of  the  tower  and  dam  and  to 
establish  some  patterns  in  the  behavior.  The  configurations  of  the  tower  and 
earthen  dam  (described  in  Chapter  3)  were  the  same  in  all  models.  The  results  of 
the  analyses  will  be  presented  in  Chapter  6  where  a  detailed  discussion  of  the 
most  important  findings  of  this  research  is  also  included. 


32 


Chapter  4  Finite  Element  Models 


5  The  MATLAB-Based 
Program 


The  primary  goal  of  this  research  was  to  perform  an  analysis  of  a  typical 
intake-outlet  tower  considering  the  multiple-support  excitation  effects.  As  men¬ 
tioned,  other  effects  such  as  the  soil-structure  interaction  were  studied  to  obtain  a 
better  understanding  of  the  dynamic  behavior  of  these  towers.  Several  tools  were 
used  to  analyze  the  towers  considering  some  or  all  of  these  effects.  One  of  the 
tools  used  during  this  research  was  a  MATLAB-based  program. 

The  finite  element  model  of  the  tower,  bridge,  dam,  and  soil  foundation 
accounts  for  the  two  effects  to  be  studied.  However,  this  model  is  relatively 
complicated  to  develop,  and  thus  a  simpler  model  was  sought.  This  model  can  be 
especially  useful  when  one  wants  to  perform  the  seismic  evaluation  of  many 
towers  in  an  inventory.  Therefore,  it  was  decided  that  MATLAB  would  be  the 
best  platform  to  develop  the  program  to  study  the  seismic  response  of  the  tower- 
bridge  system.  This  chapter  presents  a  description  of  the  formulation  used  in  this 
program  to  account  for  the  multiple-support  excitation,  and  the  hydrodynamic 
and  soil-structure  interaction  effects. 


Multiple-Support  Excitation 
Formulation 

Several  methods  and  procedures  have  been  developed  to  account  for  the 
multiple-support  excitation  effects  in  extended  structures  such  as  bridges,  dams, 
and  pipeline  systems.  Some  methods  include  the  effects  of  the  local  soil  condi¬ 
tions,  others  the  effects  of  traveling  waves,  etc.  (see  Chapter  2).  The  computer 
program  developed  for  this  study  in  MATLAB  includes  the  effects  of  multiple- 
support  excitations  by  means  of  time-history  analysis.  It  begins  with  the  formula¬ 
tion  of  the  stiffness  matrices  for  a  plane  frame  element.  The  plane  frame  element 
has  three  degrees  of  freedom  per  node  (the  axial  displacements  are  included). 

The  local  stiffiiess  matrix  of  each  element  is  given  by 


Chapter  5  The  MATLAB-Based  Program 


33 


0 

0 

A 

0 

0 

121 

61 

0 

121 

61 

61 

L 

61 

L 

~ 

41 

0 

~T 

21 

0 

0 

A 

0 

0 

121 

61 

0 

121 

61 

61 

L 

61 

L 

~ 

21 

0 

~T 

41 

where 


(3) 


E  =  modulus  of  elasticity 
L  =  length  of  the  element 
A  =  cross-sectional  area 
I  =  moment  of  inertia 


Hie  local  stifihess  matrices  are  assembled  into  die  global  stiffness  matrix 
using  tjpical  matrix  stractuial  analysis  procedures.  Hie  local  element  stiflfiiess 
matrix  can  be  transformed  into  the  global  coordinate  system  using 


K  =  T^KT 


(4) 


in  which  T  is  the  transfoimation  matrix  and  is  its  transpose  and  is  obtained 

from  followng  the  expression; 


cos  0 
-  sin  6 


T  = 


0 

0 

0 

0 


sin  9 
cos  6 
0 
0 
0 
0 


0  0  0  0 

0  0  0  0 

10  0  0 
0  cos  6  sin  6  0 

0  -sin  6  cos  9  0 

0  0  0  1 


(5) 


where  0  is  the  angle  between  the  global  horizontal  axis  and  a  local  axis  along 
the  element. 

Hie  program  also  makes  use  of  the  consistent  mass  matrix  of  plane  frame 
element,  i.e.. 
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The  consistent  matrix  includes  the  dynamic  coupling  that  exists  between  the 
transverse  displacements  and  rotations.  The  simplest  form  of  the  mass  matrix 
obtained  from  the  concentration  of  masses  on  each  node  of  the  element  (lumped 
mass  matrix)  ignores  this  effect.  This  matrix  contains  off-diagonal  terms, 
whereas  the  lumped  mass  matrix  has  only  diagonal  terms  (Paz  1997).  It  is  called 
“consistent”  because  the  same  displacement  shape  functions  that  are  used  for 
deriving  the  stifBiess  matrix  are  used  for  the  derivation  of  the  mass  matrix. 
Finally,  the  local  mass  matrix  is  assembled  into  the  global  mass  matrix  following 
the  same  procedure  for  the  assembling  of  the  global  stifl&iess  matrix. 

The  first  step  in  the  analysis  of  a  stmcture  subjected  to  multiple-support  exci¬ 
tations  is  the  determination  of  the  pseudo-static  displacements.  By  definition,  the 
pseudo-static  displacements  are  the  sum  of  the  displacements  of  fte  stmcture 
caused  by  a  slow  independent  motion  of  each  of  the  supports  of  the  stmcture.  To 
calculate  them,  it  is  necessary  to  have  the  full  stiffiiess  matrix  of  the  system.  The 
displacements  are  then  divided  into  interior  and  boundary  degrees  of  freedom. 
The  stiffiiess  and  mass  matrices  are  next  divided  in  the  same  way.  In  the  case  of 
multiple-support  excitations,  the  boundary  dofs  are  those  whose  motion  will  be 
that  of  the  ground  caused  by  the  earthquake.  For  the  tower  plus  bridge  model,  the 
program  reads  earthquake  acceleration  histories  in  the  horizontal  and  vertical 
directions  at  two  nodes.  The  remaining  dofs  are  referred  to  as  interior  dofs,  which 
in  this  case  are  also  the  pseudo-static  displacements.  The  divided  stiffness  matrix 
is  given  by  the  following  expression: 


(7) 


in  which,  m,  and  m*  are  the  vectors  of  interior  and  boundary  displacements, 
respectively.  It  will  be  assumed  tiiat  there  are  ni  interior  dofs  and  rib  boundary 
dofs.  In  the  submatrices  of  K,  the  subscripts  i  and  b  represent  the  terms  related  to 
interior  and  boundary  dofs. 

Solving  for  the  interior  (pseudo-static)  displacements  from  the  first  ni  equa¬ 
tions  leads  to 
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=-K7/K,^u^ 


(8) 


where  K7^  is  the  inverse  of  the  submatrix  that  contains  all  the  degree  of  fieedom 
associated  with  the  interior  displacements  (upper-left  matrix),  and  Kg,  is  the 
upper-right  matrix  that  results  after  the  division  of  the  total  stiffiiess  matrix. 
Introducing  a  ni  x  nb  matrix,  R  is  defined  as 

R  =  -Ki‘K,,  (9) 


in  which  die  pseudo-static  displacements  can  be  expressed  as 

Ui^Ruj,  (10) 

The  matrix  R  is  called  ftie  influence  matrix  because  it  describes  the  influence  of 
the  support  displacements  on  the  structural  displacements.  It  is  important  to 
emphasize  that  R  represents  a  matrix,  not  a  vector  as  in  the  case  of  identical  sup¬ 
port  motions.  In  the  latter  case,  R  becomes  the  vector  of  influence  coefficiente 
and  Ub  is  a  scalar,  i.e.,  the  ground  displacement.  In  the  case  of  multiple-support 
excitations,  each  element  of  is  a  displacement  time  fimction  induced  by  the 
earthquake.  Here,  each  column  of  matrix  R  represents  an  influence  vector  associ¬ 
ated  with  the  support  displacements  m*.  The  reactions  at  the  supports  can  be  cal¬ 
culated  using  the  relation  between  the  interior  displacements  Wf  and  the  prescribed 
displacements  «*.  First,  from  the  lower  part  of  the  system  of  Equation  7, 

Rf,  =  KbiU^  +  Kj,bUb  (11) 


Substituting  Ut  fix>m  Equation  8  into  1 1,  can  be  expressed  as 
R,=[K^-K„Ki’K^]u, 


(12) 


Usually  the  ground  motion  is  defined  in  terms  of  acceleration  time-histories. 
Therefore,  to  obtain  the  prescribed  displacements,  a  numerical  integration  has  to 
be  done.  The  method  for  numerical  integration  used  in  the  program  is  the  trape- 
^idal  rule.  The  program  reads  the  acceleration  time-histories  at  the  different  sup¬ 
port  locations  as  input  data.  Then  it  determines  the  displacement  time-histories  at 
each  support  of  the  structuml  system  by  numerical  integration. 


Once  tiiis  is  done,  the  program  calculates  the  dynamic  displacements  of  the 
structure.  To  do  that,  it  is  necessary  to  calculate  the  modal  participation  fector^. 
As  it  is  well  known,  these  fectors  give  a  measure  of  the  degree  to  which  the  n* 
mode  participates  in  the  response  of  the  stractural  system  and  are  given  by 


■p  _  ^b! 

'"“m 


/  =  l,...,w6 


(13) 
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Note  that,  opposite  to  the  case  of  uniform  base  motion,  there  are  ni  x  nb  partici¬ 
pation  factors.  That  is,  they  vary  not  only  with  the  mode  number  but  also  with  the 
particular  support.  The  quantity  L,a  is 

(14) 


where 

J  =  transpose  of  the  vibration  mode  (the  column  of  the  vibration 

mode  matrix) 

Mu  =  submatrix  of  the  mass  matrix  that  corresponds  to  the  terms 
associated  with  the  interior  dofs 

Ri  =  corresponding  column  in  the  influence  matrix  R  associated  with  the 
motion  of  the  support  / 

Variable  M„  in  Equation  13  is  the  generalized  modal  mass  and  is  expressed 
by 


Due  to  the  orthogonality  property  of  the  modes. 


0  ;k  ^  j 
Mj  ;k  =  J 


(15) 


(16) 


If  the  vibration  modes  are  normalized  with  respect  to  the  mass  matrix,  M„  is 
equal  to  one.  In  this  way,  the  program  normalizes  the  vibration  modes;  thus,  tiieie 
is  no  need  to  determine  M„. 

The  determination  of  the  displacements  that  consider  multiple-support  exci¬ 
tations  follows  a  procedure  similar  to  the  one  used  for  the  case  of  identical  sup¬ 
port  motions.  The  basic  difference  is  that  the  total  dynamic  response  of  the 
system  is  a  superposition  of  the  individual  responses  due  to  each  support  motion. 
To  determine  these  displacements,  it  is  necessary  to  solve  the  equations  of 
motion  transformed  in  modal  coordinates.  The  modal  equation  of  motion  is  then 
expressed  as  follows: 

rib 

;n  =  l,...,ni  (17) 

1=1 


where 

q„  =  modal  coordinate 
=  modal  damping  ratio 

(On  =  natural  fiequency  of  the  n*  vibration  mode 
nb  —  total  number  of  supports  in  the  system 
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Ugt  =  ground  acceleration  at  the  support  I 
ITie  right-hand  side  of  Equation  17  is  the  equivalent  modal  force,  given  by 

A'- (IS) 

1=1 

Each  equation  of  motion  is  solved  for  the  modal  coordinate  q„  using  the  Duhamel 
integral.  The  modal  coordinate  q„  represents  the  deformation  response  of  the 
mode  due  to  the  sum  of  all  the  support  excitations  %.  Once  the  solution  for  the 
modal  coordinates  q„  is  obtained,  the  physical  dynamic  displacements  can  be 
determined.  ITie  dynamic  displacements  at  die  interior  degrees  of  freedom  are 
obtained  combining  the  modal  coordinates,  that  is, 

«d(i)  =  T,^n9n(t)  (19) 

B=1 


where  N  is  the  total  number  of  modes  considered  in  the  analysis. 

Hie  total  displacement  response  of  the  system  is  the  combination  of  two 
parts.  It  is  the  sum  of  the  dynamic  displacements  and  the  pseudo-static  displace¬ 
ments  given  by  Equation  10.  Combining  the  two  parts  leads  to 

(20) 

In  this  equation,  %  (t)  is  die  element  of  the  vector  Ui,. 


Simplified  Procedure  for  Added 
Hydrodynamic  Masses 

It  is  recognized  that  the  seismic  response  of  structures  submerged  or  partially 
surrounded  by  water  is  influenced  by  the  surrounding  water.  Hie  interaction  with 
the  water  imposes  additional  dynamic  forces  on  the  structure  and  also  modifies 
its  djmamic  properties  (Liaw  and  Chopra  1974).  Some  of  the  effects  of  die  water 
surroimding  the  outside  and  on  the  inside  of  the  tower  include  the  reduction  of 
die  resonant  frequencies  (or  equivalendy,  the  increase  in  the  resonant  periods), 
especially  at  the  higher  modes.  Also,  in  some  case  studies,  the  hydrodynamic 
interaction  affects  the  displacements  of  the  tower;  as  a  consequence,  the  maxi¬ 
mum  shear  and  bending  moments  along  the  height  of  the  tower  also  increase,  but 
to  a  greater  extent  (Goyal  and  Chopra  1989a).  This  was  the  case,  for  example,  of 
the  Briones  Dam  intake  tower  located  east  of  San  Francisco  Bay  (Goyal  and 
Chopra  1989b). 

The  effect  of  surrounding  water  has  been  studied  for  decades  since  Wester- 
gaard  (1933)  first  proposed  the  concept  of  added  hydrod5Tiamic  masses  for  the 
seismic  analysis  of  dams  impounding  water.  The  idea  behind  this  concept  is  that 
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the  hydrodynamic  effect  of  the  fluid  is  equivalent  to  an  additional  water  pressure 
imposed  on  the  structure  due  to  a  unit  acceleration  of  the  base  of  the  structure. 

The  determination  of  the  added  mass  is  based  on  the  assumptions  that  the  struc¬ 
ture  is  rigid  and  the  water  is  incompressible  with  viscous  effects  neglected. 

It  has  been  shown  that  the  compressibility  of  the  water  is  of  no  practical  sig¬ 
nificance  in  determining  the  dynamic  response  of  slender  towers.  The  errors 
become  significant  only  for  the  higher  vibration  modes  of  squatty  towers  (Liaw 
and  Chopra  1974).  Neglecting  the  viscosity  is  also  a  common  and  accepted 
assumption  when  the  liquid  is  water.  The  effect  of  the  surface  waves  was  shown 
to  be  of  little  consequence  in  the  earthquake  response  of  towers  surrounded  by 
water  (Liaw  and  Chopra  1974).  Moreover,  the  effect  of  the  vertical  components 
of  the  ground  motion  is  ignored,  since  the  stresses  in  slender  towers  are  caused 
primarily  by  lateral  flexural  deformations. 

The  rigid  structure  assumption  was  examined  by  Liaw  and  Chopra  (1974). 
They  found  that  the  magnitude  and  distribution  of  the  additional  mass  depend  on 
the  deflected  shape  during  the  vibration.  They  developed  a  method  that  includes 
the  shape  of  the  vibration  mode  on  the  added  mass.  However,  they  had  to  con¬ 
sider  a  single  mode  to  determine  the  new  added  mass;  it  is  not  possible  to 
account  for  all  the  important  modes  of  vibration  of  the  tower.  Therefore,  they 
recommended  using  ^e  traditional  definition  of  added  mass  and,  hence,  in  subse¬ 
quent  papers,  Chopra  and  co-workers  (Goyal  and  Chopra  1989a,b,c)  adopted  this 
concept  to  study  different  aspects  of  the  seismic  response  of  intake  towers. 

With  the  assumptions  previously  mentioned,  it  is  possible  to  obtain  added 
hydrodynamic  mass  functions  that  account  for  the  water  outside  and  inside  the 
tower  for  cylindrical  towers  with  circular  cross  sections  in  a  relatively  simple 
way. 

The  expressions  for  the  added  masses  developed  by  Goyal  and  Chopra 
(1989b)  were  obtained  based  on  the  analytical  solution  of  the  Laplace  equation. 
Equation  20  shows  the  three-dimensional  Laplace  equation  in  polar  coordinates 
used  for  the  analysis  of  the  surroimding  (outside)  water  domain.  This  equation 
represents  the  irrotational  motion  of  the  water  surrounding  the  tower.  The  hydro- 
dynamic  pressure  in  the  water  surrounding  the  tower  is  generated  by  the  acceler¬ 
ation  of  the  outside  surface  of  the  tower  and  the  vertical  acceleration  of  the 
reservoir  bottom. 


d^p  ^  \  dp  ^  1  ay  I  -0  (21) 

dr^  r  dr  d6^  dz^ 


The  solution  of  the  Laplace  equation  governing  the  behavior  of  the  surround¬ 
ing  outside  water  provides  the  following  expression  for  the  corresponding  added 
mass  for  circular  cylindrical  towers  (Goyal  and  Chopra  1989b): 


\m-l 


16H^f  {-!) 


^o/hJ  cos  (a„  z/ifj 


(22) 
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where 


2  =  distance  above  the  base  of  the  tower 
Pw  =  mass  density  of  water 
To  =  radius  of  the  outside  surfece  of  the  tower 
Ho  =  depth  of  the  surrounding  winter 
a,„  =  {2m-l)nl2 


and 


Vo 


Kohm  rolH^+K^ia^  rjH^) 


(23) 


in  which  K„  is  the  modified  Bessel  function  of  order  n  of  the  second  kind. 


For  an  infinitely  long  uniform  tower  with  the  same  circular  cross  section,  the 
added  mass  per  unit  height  is 


(24) 


Note  that  wr  °  is  equal  to  the  added  m^s  of  the  water  displaced  by  the  tower  (per 
unit  height). 

Figure  15  displays  the  normalized  added  mass  ^sociated  with  the  outeide 
w^ater  as  a  function  of  for  a  range  of  values  of  jH  ^  ,  the  ratio  of  the  out¬ 

side  radius  to  the  depth  of  the  outside  water.  The  different  cinves  in  Figure  15  are 
obtained  fi’om  Equation  22  forgiven  values  of  H^jr^  and  dividing  the  results 

by .  When  the  ratio  jr^  tends  to  infinity  (diat  is,  for  an  infinitely  slender 
cjdindrical  tower),  the  normalized  added  mass  becomes  imity.  As  the  slenderness 
ratio  Hgjr^  decreases,  i.e.,  as  the  tower  becomes  squatter,  the  ratio  m°  {z)fml, 
decreases. 

Goyal  and  Chopra  (1989b)  also  studied  the  evaluation  of  added  hydrody¬ 
namic  masses  for  towers  with  noncircular  cross  sections.  Hiey  found  that  the 
normalized  ^ded  mass  for  a  uniform  tower  of  arbitrary  cross  section  is  essen¬ 
tially  the  same  as  that  of  an  equivalent  elliptical  tower.  They  defined  a  slender¬ 
ness  ratio  la^  and  a  plan  dimensions  ratio  of  the  equivalent  elliptical 

tower  in  terms  of  the  corresponding  ratios  for  fixe  actual  tower,  as  follows: 
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Figure  15.  Plot  for  normalized  outside  hydrodynamic  added  mass  for  circular  cylindrical  towers 


K 


(26) 


where 

Ho  =  depth  of  surrounding  water  for  both  towers 

2a^  =  width  of  the  elliptical  tower  perpendicular  to  the  direction  of 

ground  motion 

2b^  =  width  of  the  actual  tower  parallel  to  the  direction  of  ground 
motion 

2a^  =  width  of  the  actual  tower  perpendicular  to  the  direction  of  ground 

motion 

2b^  =  width  of  the  elliptical  tower  parallel  to  the  direction  of  ground 
motion 

Thus,  if  the  quantity  ml  I  ml^  were  available  for  towers  of  elliptical  cross 
section  for  a  range  of  values  of  the  ratio  /  b^  and  slenderness  ratio  the 

normalized  added  mass  could  be  readily  determined  for  towers  of  arbitrary  cross 
sections.  However,  Goyal  and  Chopra  (1989b)  judged  this  approach  to  be 
impractical  because  it  would  require  a  large  number  of  graphs  of  normalized 
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Figure  16.  Properties  of  “equivalent”  circular  cylindrical 
towers  corresponding  to  uniform  elliptical 
towers  associated  with  added  hydrodynamic 
mass  due  to  surrounding  water  (from  Chopra 
and  Goyal  1989a) 


added  mass  to  cover  all  the  possible  practical  values  of  I and  They 

instead  proposed  to  replace  the  elliptical  tower  by  an  equivalent  circular  tower. 
To  this  end,  they  developed  a  graph  that  gives  the  ratio  IH^  as  a  function  of 

the  ratio  of  elliptical  towers  for  different  fixed  values  of  .  This 

graph  is  displayed  in  Figure  16.  Once  an  equivalent  circular  tower  is  defined  in 
this  way.  Equation  22  or  the  graph  in  Figure  15  can  be  used  to  obtain  the  nonnal- 
ized  added  mass  of  the  equivalent  circular  tower.  Thus,  the  approach  proposed  by 
Goyal  and  Chopra  requires  defining  two  equivalent  towers:  first,  one  must  define 
a  tower  with  an  elliptical  cross  section  and,  next,  another  with  a  circular  cross 
section. 


In  summary,  the  procedure  for  defining  the  added  hydrodynamic  mass 
associated  with  the  surrounding  water  for  uniform  towers  of  arbitrary  cross  sec¬ 
tion  (with  two  axes  of  symmetry)  consists  of  five  steps: 

a.  The  parameters  ,  bo,  and  Ho  of  the  actual  tower  are  used  in  Equa¬ 
tions  25  and  26  to  determine  the  ratios  //„  /«„  and  I  b^  of  the 
equivalent  elliptical  tower.  For  a  square  cross  section,  these  equations 
yield 
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(27) 


r 


b.  Knowing  Ib^  and  the  inverse  of  the  ratio  H^la^  determined  in  Step  1, 
Figure  16  provides  the  ratio  I of  an  equivalent  circular  tower, 
where  is  its  equivalent  radius. 

c.  With  the  ratio  r^l Figure  15  allows  us  to  evaluate  the  normalized 
added  mass  ml  {z)lml^  for  the  equivalent  circular  tower.  Alternatively, 
the  ratio  ml  {z)l  ml^  can  be  directly  defined  from  Equation  22  by  setting 
the  (  p^Tirl )  equal  to  1 . 

d.  The  added  hydrodynamic  mass  for  an  infinitely  long  tower  with  the  arbi¬ 
trary  cross  section  ml^  needs  to  be  determined  using  a  table  presented  by 
Goyal  and  Chopra  (1989b-Table  1).  This  table  contains  l{p^A^)  for 
a  few  selected  cross  sections.  Based  on  the  data  in  this  table,  a  graph 
depicting  the  relationship  between  the  ratio  Qolbo  and  the  normalized 
added  hydrodynamic  mass  ml^l{p^A^)  for  a  rectangular  section  is  pre¬ 
sented  in  Figure  17.  For  a  square  cross  section  {qo  =  b^,  the  ratio 

</(pA)  is  1.186. 


Figure  17.  Normalized  outside  hydrodynamic  added  mass  for  an  infinitely 
long  tower 
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e.  Finally,  the  added  hydrodynamic  mass  ml  (z)  for  the  actual  tower  is 
computed  by  multiplying  the  normalized  added  mass  ml  iz)lml^  calcu¬ 
lated  in  Step  3  by  w®  defined  in  Step  4. 


Similarly  to  the  expression  for  the  case  of  the  surrounding  water,  it  can  be 
shown  that,  firom  the  solution  of  the  Laplace  equation  for  the  inside  water  for 
circular  cylindrical  towers,  one  obtains 


W=( 


^AELfMl 


COS  la^Zi 


where 


z  =  distance  above  the  base  of  the  tower 
Vi  -  radius  of  the  inside  surface  of  the  tower,  =  {2m  -  i)^/ 2 

Hi  =  depth  of  the  inside  water 


I O («,«  rt iHi )+ 1 2{a„rJ Hi ) 


(29) 


in  which  is  the  modified  Bessel  function  of  order  n  of  the  first  kind. 


The  normalized  added  mass  is  obtained  fi-om  Equation  28  dividing  wa^(z)  by 
m'  -  p  7tr^ ,  Figure  18  displays  the  normalized  added  mass  associated  with  the 

water  contained  inside  the  tower  for  a  range  of  values  of  jH^ ,  the  ratio  of  the 
inside  radius  to  inside  water  depth.  When  the  ratio  tends  to  infinity  (i,e., 

for  an  infinitely  slender  cylindrical  tower),  the  normalized  added  mass  is  unity 
and  it  decreases  as  the  tower  becomes  squatter  (i.e,,  the  ratio  decreases). 

Comparing  the  normalized  added  mass  for  outside  water  with  the  inside  water 
masses,  it  is  found  that,  for  the  same  slenderness  ratio,  the  normalized  inside 
added  mass  is  larger. 

For  a  circular  cylindrical  tower,  Equations  28  and  29  completely  define  the 
added  hydrodynamic  mass.  The  equations,  however,  are  not  directly  applicable 
for  towers  with  noncircular  cross  sections.  Goyal  and  Chopra  (1989b)  also  devel¬ 
oped  a  simplified  procedure  to  calculate  the  added  hydrodynamic  mass  for 
towers  with  arbitrary  cross  section,  provided  they  have  two  axes  of  symmetry. 

Let  the  2a/  and  21?/  be,  respectively,  the  interior  dimensions  of  the  actual  tower 
perpendicular  and  along  the  direction  of  the  ground  motion  Hj  (the  interior  water 
depth),  and  ^/,  the  interior  cross-section  area  of  the  tower.  The  methodology  is 
similar  to  the  one  presented  before  and  consists  of  the  following  steps: 
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Figure  18.  Plot  for  normalized  inside  hydrodynamic  added  mass  for  circular 
cylindrical  towers 


a.  Obtain  the  dimensions  of  an  equivalent  elliptical  tower  with  interior 
cross-sectional  dimensions  2a.  and  2b.  in  terms  of  the  dimensions  of  the 
actual  tower,  as  follows: 


(30) 


(31) 


For  a  tower  with  a  square  cross  section,  these  equations  yield 


H., 


-n 


a. 


1 


(32) 


h.  Calculate  the  slenderness  ratio  H.  /  r.  of  an  equivalent  circular  tower 
using  the  properties  determined  in  Step  1.  According  to  Goyal  and 
Chopra  (1989b),  for  the  inside  water  analysis,  the  equivalent  radius  r  is 
given  by  the  following  simple  formula,  instead  of  a  graph  as  in  the  case 
of  the  surrounding  water: 
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(33) 


and  r.  =  for  a  tower  with  a  square  cross  section. 

c.  With  H.  and  7}  ,  evaluate  the  normalized  added  mass  for  the 

equivalent  cylindrical  tower  using  Figure  18.  In  this  case,  the  added 
hydrodynamic  mass  for  the  tower  with  arbitrary  cross  section  is  obtained 
by  multiplying  the  normalized  added  mass  by  , 

d  As  an  alternative  to  Step  c,  the  added  hydrodynamic  mass  w*  (z)  can  be 
obtained  directly  from  Equation  28. 

Note  that  in  the  case  of  the  inside  water,  is  the  same  regardless  of  the 
shape  of  the  tower's  cross  section.  In  other  words,  it  is  not  necessary  to  have 
available  tables  or  graphs  to  determine  the  added  hydrodynamic  mass  for  an 
infinitely  high  tower  with  different  cross  sections,  as  in  the  case  of  the  surround¬ 
ing  water. 

Once  the  added  hydrodynamic  masses  for  the  outside  m^(z)  and  inside 
water  are  obtained,  the  next  step  is  to  lump  these  distributed  masses  and 
add  them  to  the  discrete  masses  due  to  the  inertia  of  the  tower.  The  procedure 
varies  depending  on  whether  a  shell  or  frame  model  is  used  to  represent  the 
tower.  In  the  case  of  a  frame  model  such  as  the  one  used  in  this  study,  the  proce¬ 
dure  is  the  following.  A  “tributary’'  length  must  be  defined  for  each  lumped  mass 
along  the  tower  height,  say  Lj  for  the  mass  at  the node.  Suppose  that  the 
node  is  at  a  height  hi  measured  from  the  base  of  the  tower.  Then,  the  lumped 
mass  at  the  node  due  to  the  outside  and  inside  water  is  approximately  given  by 

Mj=[m:ihj)  +  m:ihj)]Lj  (34) 

A  more  accurate  value  would  be  obtained  by  using 


However,  this  requires  performing  a  very  complicated  integration  involving  spe¬ 
cial  functions.  Therefore,  Equation  34  will  be  used  in  die  present  study. 


Soil-Structure  Interaction 

Another  factor  that  has  a  strong  influence  on  the  seismic  response  of  intake- 
outlet  towers  is  the  soil-interaction  effect.  Initially  this  effect  was  not  going  to  be 
considered.  It  is  well  known  that  the  soil-structure  interaction  effects  are 
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important  only  for  stiff  and  massive  structures  located  on  soft  soils,  and  it  was 
thought  the  intake  tower  did  not  comply  with  these  characteristics.  However,  the 
discrepancies  between  the  simplified  frame  model  and  the  fiill  finite  element 
model  of  the  tower  plus  dam  (reported  in  Chapter  4)  led  to  its  inclusion. 

The  most  common  approach  to  the  problem  is  to  modify  the  dynamic  proper¬ 
ties  of  the  structure  and  evaluate  the  response  of  the  modified  structure  to  the 
prescribed  ffee-field  ground  motion  (Building  Seismic  Safety  Council,  BSSC 
2001).  Among  the  approaches  to  achieve  this,  the  simplest  one  involves  consider¬ 
ing  the  flexibility  of  the  soil  by  means  of  translational  and  rotational  springs. 
Moreover,  because  a  substantial  part  of  the  structure’s  vibrational  energy  is  dissi¬ 
pated  into  the  supporting  medium  by  radiation  of  energy  and  hysteretic  behavior 
of  the  soil,  translational  and  rotation  dashpots  are  added  to  the  model.  The  first 
effect,  called  radiation  damping,  is  typically  much  more  important  than  the  sec¬ 
ond,  which  is  referred  to  as  material  damping. 

The  method  for  incorporating  the  effects  of  soil-structure  interaction  based 
on  springs  and  dashpots  was  thought  to  provide  sufficient  accuracy  for  the  appli¬ 
cation  at  hand  and,  thus,  was  the  approach  adopted  in  this  study.  The  coefficients 
of  the  springs  and  dashpots  were  calculated  using  the  expressions  described  in 
Chapter  4  or  the  values  obtained  by  using  SAP2000  (explained  in  Chapter  6).  As 
recommended  by  Newmark  and  Rosenblueth  (1971),  lumped  masses  were  also 
attached  to  the  degrees  of  freedom  at  the  soil-structure  interface.  More  elaborate 
procedures  have  their  own  approximations  and  do  not  eliminate  the  uncertainties 
that  are  inherent  in  the  modeling  of  the  structure-foundation-soil  system,  in  the 
specifications  of  the  ground  motion,  and  in  the  estimation  of  the  properties  of  the 
stmcture  and,  especially,  the  soil  (BSSC  2001). 

The  effect  of  soil-structure  interaction  at  the  base  of  the  tower  and  at  the 
location  of  the  bridge  was  introduced  in  the  MATLAB  program  modifying  the 
stiffness  and  mass  matrices  of  the  structural  system.  The  values  of  the  mass  and 
stifrhess  coefficients  defined  in  Chapter  4  were  added  to  the  degrees  of  freedom 
at  the  interface  between  the  structure  and  the  soil  foimdation.  Because  the  springs 
representing  the  flexibility  of  the  soil  have  one  end  attached  to  a  fixed  point,  their 
stif&iess  coefficients  are  added  to  die  diagonal  terms  of  the  stiffness  matrix  at  the 
base  of  the  tower  and  at  one  end  of  the  bridge.  The  same  occurs  with  the  mass 
coefficients. 

The  radiation  and  material  damping  of  the  soils  are  accounted  for  by  means 
of  discrete  dashpots.  As  explained  in  Chapter  4,  these  dashpots  were  defined  only 
for  the  base  of  the  tower  because  in  the  technical  literature  they  are  determined 
only  for  the  case  of  an  elastic  half-space.  In  principle,  the  coefficients  of  the 
dashpots  in  the  three  directions  (horizontal,  vertical,  and  rotational  or  rocking) 
should  be  added  to  the  damping  matrix  of  the  structure.  However,  there  is  a 
problem  in  this  case:  the  damping  matrix  of  the  structure  was  not  defined.  Rather, 
the  usual  classical  damping  assumption  was  used,  and  the  damping  of  the  struc¬ 
ture  was  introduced  by  means  of  modal  damping  ratios.  In  other  words,  if  O  is 
the  matrix  of  eigenvectors  of  the  tower  (or  tower  plus  bridge,  depending  on  the 
model  considered)  with  a  flexible  foundation  (i.e.,  with  the  springs  at  the  base),  it 
is  assumed  that  the  triple  product 
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0^C0  =  fftag[c]  = 


(36) 


leads  to  a  diagonal  matrix.  Furthermore,  for  convenience,  the  diagonal  terms  are 
expressed  in  terms  of  the  naturd  frequencies  ©j  and  modal  (temping  ratios  §j: 

c.='24p.  ;  J  =  (37) 

If  no  other  information  is  available,  the  values  of  the  damping  ratios  §j  are  usu¬ 
ally  considered  to  be  equd  to  0.05  in  order  to  use  design  spectra  specified  in 
codes  or  seismic  provisions.  However,  for  consistency  with  the  finite  element 
model  of  the  dam  plus  inlake  tower,  the  value  of  §  was  set  equal  to  0.245.  It  is 
recalled  that  this  value  accounts  for  the  nonlinear  behavior  of  the  soil  deposit  and 
was  derived  by  matching  die  surfece  ^celeratioiK  calculated  with  QUAD4M  and 
SAP2000. 

When  the  dashpots  are  added  to  the  base  of  the  tower,  the  total  damping 
matrix  Q  is  now  the  sum  of  two  terms: 

C,=C+Q  (38) 

where  Cj  is  a  damping  matrix  with  only  tiiree  nonzero  terms  in  correspondence 
wth  the  two  translational  and  rotational  dofe  at  the  base  of  the  tower.  When  the 
new  equations  of  motion  are  decoupled  using  the  eigenvector  matrix  the  usual 
way,  die  product 

=  (39) 

does  not  lead  to  a  diagonal  matrix.  That  is,  is,  in  general,  a  full  matrix.  In  prin¬ 
ciple,  this  situation  can  be  taken  into  account  in  an  exact  way  by  transforming  the 
usual  equations  of  motion  into  “state  space”  equations  and  solving  them  using 
complex  eigenvalues  and  eigenvectors.  However,  the  implicit  imceitainties  in  the 
values  of  the  damping  coefficients  and  the  assumptions  and  simplifications 
adopted  so  fer  do  not  warrant  the  (implications  introduced  in  the  analysis  by  this 
procedure. 

Taking  into  acciimt  that  the  diagonal  terms  of  are  dominant  with  respect 
to  the  off-diagonal  terms,  it  is  proposed  to  neglecrt  the  latter  terms  and  replace 
Equation  39  with 


(40) 
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With  this  assumption,  the  total  damping  matrix  C,  can  now  be  decoupled  using 
the  modes  of  the  tower  with  flexible  base,  that  is 


=  diag]c  ] + diag\c^  ]  = 


C+C- 

j  sj 


(41) 


With  the  diagonal  terms  in  Equation  30,  one  can  define  new  equivalent  modal 
damping  ratios  that  approximately  account  for  the  presence  of  the  dashpots,  as 
follows: 


And,  using  Equation  26,  one  can  further  write 


S' 


(42) 


(43) 


The  new  modal  damping  ratios  are  used  in  the  traditional  modal  analysis  to 

obtain  the  seismic  response  of  the  intake  tower  with  soil-stmcture  interaction 
effects  included. 


Summary 

A  description  of  the  formulation  used  in  the  MATLAB  program  was  pre¬ 
sented  in  this  chapter.  The  following  features  are  included  in  the  program:  the 
multiple-support  excitation  formulation  and  simplified  procedures  to  define  the 
added  hydrodynamic  masses  and  the  soil-structure  interaction  effects.  The  prog¬ 
ram  receives  all  the  data  related  to  the  structure  in  an  interactive  way.  The  output 
of  the  program  is  the  absolute  acceleration  and  relative  displacement  time-history 
of  any  mode  of  interest. 

The  first  step  in  considering  the  multiple-support  excitations  is  the  calcula¬ 
tion  of  the  pseudo-static  displacements.  These  displacements  are  obtained  by 
applying,  in  a  pseudo-static  way,  prescribed  displacements  at  each  support  of  the 
stmctural  system.  In  the  case  of  multiple-support  excitations,  the  prescribed  dis¬ 
placements  are  those  induced  by  the  earthquake.  Afterwards,  the  dynamic  dis¬ 
placements  are  calculated  using  the  usual  modal  analysis  formulation.  The  total 
response  is  the  sum  of  the  pseudo-static  and  dynamic  displacements. 

Another  feature  included  in  the  program  is  the  simplified  procedure  to  define 
the  added  hydrodynamic  masses  developed  by  Goyal  and  Chopra  (1989b).  This 
procedure  is  used  to  obtain  the  Ivunped  masses  that  accoimt  for  the  presence  of 
the  water  inside  and  outside  the  tower.  The  hydrodynamic  added  masses  were 
added  to  the  masses  of  the  tower  of  the  frame  model  and  in  the  tower  of  all  the 
finite  element  models. 
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The  last  feature  included  in  the  program  is  the  soil-stracture  interaction 
effects.  Hiis  effect  is  included  in  a  simple  way  using  springs,  raises,  and  dash- 
pots  that  represent  the  soil  foundation  beneath  the  tower  and  in  the  dam  w^here 
the  bridge  is  supported.  The  springs  and  masses  are  simply  accounted  for  by 
modifjdng  the  appropriate  terms  of  the  stififiiess  and  mass  matrices.  The  Edition 
of  dashpots  requires  an  additional  assumption  in  order  to  avoid  the  use  of  com¬ 
plex  eigenvalues  and  eigenvectors.  Hie  off-diagonal  terms  in  the  matrix  obtained 
by  pro-  and  post-multiplying  the  damping  matrix  containing  the  dashpols  are 
neglected,  and  new  modal  damping  ratios  are  calculated  with  diagonal  terms. 
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6  Numerical  Results 


This  chapter  presents  a  summary  of  the  results  obtained  during  the  investiga¬ 
tion.  The  results  are  divided  in  two  groups:  those  obtained  with  the  finite  element 
models  and  those  obtained  with  the  frame  model  developed  in  MATLAB.  The 
results  of  the  seismic  analysis  are  displayed  in  graphical  form.  Only  the  relative 
displacement  and  absolute  acceleration  time-histories  are  presented.  To  reduce 
the  number  of  plots  presented  in  the  report,  only  the  response  at  the  top  of  the 
tower  is  shown.  However,  the  data  are  available  for  different  points  of  the  struc¬ 
tural  system.  The  results  of  an  eigenvector  analysis  are  presented  in  tables.  A 
comparison  of  responses  obtained  with  the  different  models  and  assumptions  is 
also  presented.  Different  results  are  obtained  by  either  including  or  neglecting  the 
factors  that  can  affect  the  seismic  response  of  the  intake-outlet  towers.  These 
factors  include  the  multiple-support  excitation  and  the  soil-structure  interaction 
effects.  A  detailed  discussion  of  the  most  important  findings  of  this  research  is 
also  presented. 


Seismic  Input 

The  first  step  in  any  seismic  analysis  is  the  definition  of  the  seismic  input.  In 
this  study,  the  seismic  input  for  all  the  finite  element  and  frame  models  was 
defined  in  terms  of  an  accelerogram  applied  at  the  base  of  the  dam’s  foundation 
(at  the  bedrock).  A  response  spectrum  analysis  was  also  carried  out,  and  the 
results  are  presented  at  the  end  of  the  chapter.  The  results  are  first  presented  for 
two  earthquakes.  (Later  on,  more  ground  motions  are  considered.)  The  first  one 
is  the  classical  El  Centro  earthquake  recorded  in  Imperial  Valley,  California,  on 
15  October  1940.  Figure  19  shows  a  plot  of  the  acceleration  time-history  of  the 
El  Centro  earthquake  used  in  this  chapter.  All  the  accelerations  are  presented  as  a 
fraction  of  g  (the  gravitational  acceleration).  For  brevity,  this  is  indicated  in  the 
vertical  axes  as  “%g.”  The  peak  acceleration  of  the  El  Centro  earthquake  is 
0.319  g,  and  it  has  a  total  duration  of  30  sec.  This  accelerogram  was  used  because 
it  became  a  standard  input  in  most  earthquake  engineering  studies  and  also 
because  it  represents  a  typical  strong,  broad-band  motion.  The  other  ground 
motion  selected  for  the  analyses  is  the  Loma  Prieta  earthquake,  recorded  at  the 
Gilroy  Historic  Building,  in  California,  on  18  October  1989.  Figure  20  presents 
its  acceleration  time-history,  which  has  a  total  duration  of  40  sec  with  a  peak 
acceleration  of  0.284  g.  This  accelerogram  was  used  because  it  has  a  different 
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Figure  22.  Smoothed  Fourier  spectrum  of  the  Loma  Prieta  earthquake  at 
Gilroy  Building 
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Figure  23.  Raw  and  smoothed  Fourier  spectra  for  El  Centro  earthquake 


Figure  24.  Raw  and  smoothed  Fourier  spectra  for  Loma  Prieta  earthquake  at 
Gilroy  Building 
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Figure  25.  Arias  intensity  for  El  Centro  earthquake 


Time  (sec) 


Figure  26.  Arias  intensity  for  Loma  Prieta  earthquake  at  Gilroy  Building 
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Figures  23  and  24  display  the  raw  and  smoothed  Fourier  spectra  in  terms  of 
frequency  for  the  two  accelerograms  selected  for  this  study.  The  original  Fourier 
spectra  of  both  earthquakes  are  very  jagged,  but  when  the  smoothed  spectra  are 
plotted,  the  differences  become  apparent.  Note  that  the  Fourier  spectrum  for 
Loma  Prieta  has  more  pronounced  peaks  at  low  frequencies. 

Figures  25  and  26  show  the  Arias  intensity  curves  used  for  the  calculation  of 
the  duration  of  the  strong  motion  of  each  earthquake  (El  Centro  and  Loma 
Prieta).  It  can  be  seen  that  the  El  Centro  earthquake  has  higher  longer  motion 
duration  than  the  Loma  Prieta  earthquake. 


Multiple-Support  Excitation 
Analysis 

The  multiple-support  excitation  effects  in  the  intake-outlet  tower  and  access 
bridge  system  were  studied  using  the  tools  described  in  the  previous  chapters. 

The  finite  element  model  of  the  dam  plus  the  intake  tower  and  the  frame  model 
developed  in  MATLAB  (using  the  formulation  described  in  Chapter  5)  were  used 
in  this  phase  of  the  study.  The  finite  element  model  of  the  dam  without  the  tower 
was  used  to  obtain  the  absolute  acceleration  at  the  base  of  tower  and  at  the  point 
at  which  the  access  bridge  is  supported  by  the  dam.  The  support  of  the  bridge 
(i.e.,  the  abutment)  is  assumed  to  be  at  the  crest  of  the  embankment.  These  accel¬ 
eration  time-histories  were  used  as  input  for  the  MATLAB  program.  Figures  27 
and  28  show  a  comparison  between  the  acceleration  time-histories  at  the  base  of 
the  tower  and  at  the  access  bridge  for  the  El  Centro  and  Loma  Prieta  earthquakes, 
respectively.  These  accelerograms  were  obtained  from  the  finite  element  model 
of  the  dam  plus  foundation  created  in  SAP2000. 

The  variability  in  the  acceleration  at  the  base  of  the  tower  and  at  the  crest  of 
the  dam  (at  the  abutment  of  the  access  bridge)  is  less  than  was  expected,  at  the 
beginning  of  the  study.  There  is,  though,  a  slight  amplification  at  the  top  (crest) 
of  the  dam.  The  lack  of  a  more  pronounced  difference  could  be  the  result  of  the 
geometry  and  soil  properties  of  the  dam.  The  ratio  between  the  dam^s  height  and 
the  width  of  its  base  is  relatively  small  (0.14).  The  materials  of  the  dam  have 
shear  wave  velocities  ranging  from  450  to  700  ft/sec. 

The  two  accelerograms  for  each  earthquake  (shown  in  Figures  27  and  28) 
were  used  as  input  to  the  tower  plus  bridge  frame  model.  The  acceleration  at  the 
base  of  the  tower  was  then  applied  to  a  stick  model,  and  the  relative 
displacements  at  the  top  of  the  tower  were  compared  for  the  two  cases. 

Figures  29  and  30  display  the  relative  displacement  time-histories  at  the  top 
of  the  tower  obtained  by  considering  the  multiple-support  excitations  (MSE)  and 
a  single  excitation  (SE).  The  two  sets  of  results  were  obtained  with  the  frame  and 
column  model  of  the  intake  tower  implemented  in  MATLAB  and  using  the  for¬ 
mulation  discussed  in  Chapter  5  for  the  MSE  case.  The  results  demonstrate  that, 
at  least  for  the  case  studied,  the  effect  of  multiple-support  excitations  is  not  sig¬ 
nificant  because  both  displacements  are  very  similar.  The  amplitudes  of  the 
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Figure  27.  Acceleration  at  the  tower  base  and  bridge  support  (El  Centro) 


Figure  28.  Acceleration  at  the  tower  base  and  bridge  support  (Loma  Prieta) 
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Figure  30.  Relative  displacements  (Loma  Prieta) 
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displacements  calculated  with  the  two  models  have  a  small  variation  along  the 
complete  time-history.  Only  a  small  difference  in  the  frequency  content  of  the 
response  can  be  noted.  Thus,  at  least  for  this  particular  case,  the  effect  of  the 
spatial  variability  of  the  ground  motion  at  the  supports  is  not  important.  This 
result  was  expected  because  of  the  small  variation  in  the  acceleration  time- 
histories  at  each  support.  It  is  possible,  of  course,  that  in  other  situations  in  which 
the  seismic  motion  at  the  top  of  the  dam  is  significantly  amplified  and  filtered, 
the  consideration  of  the  MSE  effects  could  be  important. 

To  gain  further  insight  into  the  response  of  the  structural  system,  the 
response  was  also  calculated  with  the  finite  element  model  of  the  dam  plus  the 
intake  tower  and  bridge.  This  permits  verifying  the  results  of  the  MSE  analysis 
with  a  more  accurate  (but  more  complicated)  model.  Figures  31  and  32  show  a 
comparison  between  the  relative  displacement  time-histories  at  the  top  of  the 
tower  predicted  by  the  dam  plus  tower  finite  element  model  and  by  the 
MATLAB-based  frame  model  with  the  MSE. 
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Figure  31 .  Relative  displacements  obtained  with  the  frame  model  with  MSE  effects  and 
with  the  dam  plus  tower  finite  element  model  (El  Centro) 


Obviously,  the  results  are  extremely  different,  almost  of  one  order  of  magni¬ 
tude.  At  first  sight,  these  results  were  surprising  since  it  was  previously  shown 
that  the  MSE  effects  were  not  important.  These  findings  demonstrated  the  neces¬ 
sity  of  carrying  out  more  analyses  and  extending  the  investigation. 

To  begin  the  investigation  into  the  causes  of  the  differences,  the  acceleration 
time-histories  at  the  base  of  the  tower  and  at  the  support  of  the  bridge  were  now 
calculated  using  the  finite  element  model  of  the  dam  with  the  intake  tower 
included.  It  is  recalled  that  the  accelerations  used  as  input  to  the  frame  model 
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Figure  32.  Relative  displacements  obtained  with  the  frame  model  with  MSE  effects  and 
with  the  dam  plus  tower  finite  element  model  (Loma  Prieta) 


were  obtained  from  the  finite  element  model  of  the  dam  alone  (presented  in  Fig¬ 
ures  27  and  28).  The  absolute  accelerations  at  the  base  of  the  tower  and  at  the 
bridge  abutment  (at  the  crest  of  the  dam)  computed  with  the  combined  model  are 
displayed  in  Figures  33  and  34  for  the  El  Centro  and  Loma  Prieta  earthquakes, 
respectively.  It  is  very  interesting  to  note  the  difference  between  the  acceleration 
time-histories  predicted  by  these  two  models  (with  and  without  the  tower).  The 
first  difference  that  can  be  noticed  is  that  the  magnitudes  of  the  accelerations  are 
quite  different.  The  accelerations  will  be  compared  in  the  following  section, 
where  additional  discussion  will  also  be  presented.  Here  we  would  like  to  point 
out  an  interesting  and  unexpected  phenomenon.  The  majority  of  the  previous 
investigations  dealing  with  dam  response  reported  an  amplification  of  the  accel¬ 
erations  at  the  top.  However,  the  results  presented  in  Figures  33  and  34  reveal 
that  the  situation  is  different  when  the  intake  tower  and  the  bridge  are  included  in 
the  model  of  the  dam  and  its  foundation.  Here,  there  is  a  slight  attenuation  of  the 
acceleration  amplitudes  at  the  crest  of  the  dam.  It  can  also  be  observed  from  the 
figures  that  the  acceleration  time-histories  at  the  top  of  the  dam  (at  bridge  sup¬ 
port)  have  the  higher  frequency  components  filtered  out,  as  expected. 


Soil-Structure  Interaction 

To  provide  a  better  picture  of  the  differences  between  the  accelerations  act¬ 
ing  on  the  tower  base  obtained  with  finite  element  models  of  the  dam  alone  and 
the  dam  plus  tower,  two  comparative  plots  are  presented.  Figures  35  and  36  show 
the  absolute  acceleration  time-histories  at  the  base  of  the  dam  for  the  two  finite 
element  models  and  for  the  two  earthquakes  being  considered.  These  plots 
indicate  that  there  are  significant  differences  between  the  accelerations  at  the 
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Figure  33.  Accelerations  at  the  tower  base  and  bridge  support  obtained  with  the  dam  plus 
intake  tower  finite  element  model  (El  Centro) 


Figure  34.  Accelerations  at  the  tower  base  and  bridge  support  obtained  with  the  dam 
plus  intake  tower  finite  element  model  (Loma  Prieta) 


dam  base  when  the  tower  is  present  and  when  it  is  not.  The  acceleration  at  the 
dam  plus  tower  results  demonstrates  the  importance  of  the  dynamic  interaction 
between  the  intake  tower  and  the  soil  structure  (the  dam  and  the  foundation). 
These  findings  explain  the  great  differences  in  the  previous  results. 
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Figure  36.  Accelerations  at  the  base  of  the  dam  obtained  with  the  dam  plus  intake  tower 
and  the  dam  alone  finite  element  models  (Loma  Prieta) 


Another  result  obtained  with  the  model  of  the  dam  with  the  tower  that  can 
help  to  gain  insight  into  the  phenomenon  is  the  variation  of  the  free^field  acceler¬ 
ation  at  different  points  along  the  free  surface.  Three  points  were  selected  to 
compute  the  acceleration;  diese  are  shown  in  Figure  37.  Figures  38  and  39  show 
the  acceleration  time-histories  at  the  three  points  in  the  free  surface.  The  acceler¬ 
ation  time-histories  at  nodes  1 172  and  1055  shown  in  Figure  38  look  very  simi¬ 
lar.  These  nodes  are  at  the  same  elevation,  but  on  opposite  sides  of  the  dam.  This 
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Figure  38.  Acceleration  time-histories  at  two  nodes  at  the  free  surface  on  opposite  sides  of 
the  dam 


similarity  was  expected  since  the  free-field  accelerations  at  a  point  with  the  same 
elevation  on  a  soil  deposit  with  similar  conditions  should  be  very  similar.  How¬ 
ever,  as  shown  in  Figure  39,  the  acceleration  trace  at  node  1  looks  very  different 
from  that  at  node  1172,  which  is  located  just  to  the  side  of  it  (see  Figure  37). 
Node  1  is  located  at  the  interface  of  the  dam  and  the  intake  tower.  This  result 
once  again  demonstrates  the  importance  of  the  soil-structure  interaction  effects  in 
this  particular  case  study.  The  difference  between  the  motions  of  the  soil  surface 
and  of  the  base  of  the  structure  will  result  in  a  quite  distinct  response  of  the 
intake  tower  when  its  base  is  assumed  fixed. 

An  additional  numerical  test  was  done  to  assess  the  importance  of  the  soil 
structure-interaction  effects  in  this  particular  case.  The  modulus  of  elasticity 
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Figure  39.  Acceleration  time-histories  at  a  node  at  the  tower  base  and  on  the  free  surface 


of  each  material  that  composes  the  embankment  was  increased  so  that  the  point 
of  support  of  the  tower  and  bridge  become  closer  to  a  fixed-base  condition  in 
which  the  soil-structure  interaction  effects  should  diminish.  As  expected,  as  the 
moduli  were  increased,  the  acceleration  time-histories  became  closer  to  those  of 
the  dam  alone  where  no  interaction  occurs.  The  ample  differences  in  the  stiffness 
of  the  structure  compared  with  the  soil  stiffness  are  responsible  for  the  important 
soil-structure  interaction  effects.  The  values  of  the  shear  moduli  of  the  different 
types  of  soils  in  the  embankment  are  very  small  compared  with  those  of  more 
competent  soils.  The  soils  in  the  embankment  classify  as  Sd  or  Se  according  to 
the  1997  Unifonn  Building  Code  (ICBO  1997),  where  Sb  is  rock  or  competent 
soil.  Moreover,  it  is  recalled  that  the  intake  tower  is  made  of  a  hollow  concrete 
square  section  with  outside  dimensions  of  30  x  30  ft,  wall  thickness  of  3  ft,  and 
height  of  75  ft.  Therefore,  the  tower  is  very  rigid  compared  with  the  soil  deposit 
beneath  it.  As  a  result,  the  foundation  will  translate  and  rotate,  causing  the 
tower’s  displacements  to  increase  with  respect  to  those  previously  calculated  (for 
example,  those  displayed  in  Figures  29  and  30). 

To  account  for  the  soil-structure  interaction  effects,  a  simplified  model  of  the 
intake-outlet  tower  was  proposed  in  which  the  soil  stiffness  and  damping  were 
represented  by  translational  and  rotational  springs  and  dashpots,  respectively.  As 
mentioned  in  Chapter  4,  the  stiffiiess  and  damping  coefficients  were  calculated 
from  the  expressions  presented  in  the  SAP2000  reference  manual  (CSI 2000). 
These  equations  are  similar  to  those  first  proposed  by  Lysmer  (1965)  and  other 
authors  (Gazetas  1991),  which  are  independent  of  the  frequency  of  the  excitation. 
These  formulas  were  used  only  at  the  base  of  the  tower;  i.e.,  they  represent  the 
soil  below  the  foundation  of  the  tower.  The  corresponding  stiffness  coefficients 
for  the  springs  at  the  point  where  the  bridge  is  lying  on  the  dam  were  obtained  in 
a  different  way.  They  were  calculated  using  the  approach  based  on  the  applica¬ 
tion  of  a  unit  force  in  a  finite  element  model.  Using  SAP2000,  a  unit  force  was 
applied  at  the  point  where  the  bridge  is  connected  to  the  dam.  The  force  (or 
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moment)  was  applied  once  at  a  time  in  three  directions  (horizontal,  vertical,  and 
rotation),  and  the  displacements  and  rotation  along  these  directions  were 
obtained.  The  stiffness  coefficients  were  then  calculated  with  the  simple 
expression. 


in  which  F  and  A  represent  the  unit  force  (or  moment)  and  displacement  (or  rota¬ 
tion)  at  each  direction,  respectively. 

No  dashpot  was  placed  at  the  bridge  support  since  the  formulations  available 
in  the  literature  are  based  on  a  rigid  disk  in  a  semi-infinite  half-space  and  thus 
they  are  not  valid  for  the  crest  of  the  dam.  The  derivation  of  expressions  for  the 
damping  coefficients  for  this  particular  case  is  out  of  the  scope  of  this  report. 
Thus,  it  was  assumed  that  the  springs-dashpots  but  without  dashpots  at  the  bridge 
location  could  approximate  reasonably  well  the  flexibility  and  energy  radiation  of 
the  deposit  below  the  tower  and  the  dam  below  the  bridge  support.  This  is  a  rea¬ 
sonable  assumption  because  the  available  damping  coefficients  represent  the 
energy  dissipation  by  radiation  from  the  support  into  an  elastic  half-space  with  a 
plane  interface.  The  crest  of  the  dam  where  the  bridge  is  supported  cannot  be 
considered  as  an  elastic  half-space,  and  thus  the  radiation  damping  should  be 
smaller  (albeit  not  zero,  as  it  is  being  assumed  here).  Ignoring  this  damping  will 
also  lead  to  conservative  results. 

Using  the  values  of  the  coefficients  and  the  assumptions  previously 
described,  the  frame  model  was  analyzed  again  with  the  free-field  acceleration  as 
input.  The  displacements  at  the  top  of  the  tower  were  also  calculated  with  the 
dam  plus  intake  tower  model  to  compare  the  results.  Figures  40  and  41  present 
the  relative  displacement  time-histories  at  the  top  of  the  tower  for  the  frame  with 
spring-dashpots  and  for  the  dam  plus  tower  finite  element  model  for  the 
El  Centro  and  Loma  Prieta  earthquakes. 

Figures  40  and  41  show  that  the  relative  displacements  of  the  proposed 
spring-dashpot  model  and  the  dam  plus  intake  tower  finite  element  model  are  in 
good  agreement  for  both  earthquakes.  Therefore,  based  on  these  results  one  can 
preliminarily  say  that  the  simplified  spring-dashpot  model  is  reliable  for  this 
particular  case.  However,  and  in  order  to  confirm  the  statement,  additional  tests 
are  required.  Both  models  were  then  analyzed  by  subjecting  them  to  earthquakes 
with  shorter  and  longer  durations.  Figures  42  and  43  show  the  relative  displace¬ 
ments  at  the  top  of  the  tower  for  two  short-duration  earthquakes,  namely 
Coalinga  and  Anza  ground  motions.  The  Coalinga  accelerogram  was  recorded  at 
Oil  Fields  Fire  Station,  California,  on  2  May  1983  and  the  Anza  record  was 
measured  at  Pinyon  Flat  on  30  October  2001. 

These  figures  show  that  again  there  are  significant  differences  between  the 
relative  displacements  at  the  top  of  the  tower  predicted  by  the  two  models.  These 
results  were  unexpected  because  it  was  deemed  that  the  spring-dashpot  model 
would  give  a  good  approximation  for  all  seismic  motions.  Thus,  additional  tests 
were  performed  to  explain  these  findings. 
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Figure  40.  Displacements  at  the  top  of  the  tower  for  the  dam  plus  tower  and  the  frame 
with  spring-dashpot  models  (El  Centro) 


Figure  41 .  Displacements  at  the  top  of  the  tower  for  the  dam  plus  tower  and  the  frame 
with  spring-dashpot  models  (Loma  Prieta) 
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Figure  42.  Displacements  at  the  top  of  the  tower  for  the  dam  plus  tower  and  the  frame  with 
spring-dashpot  models  (Coalinga) 


— ^  Dam  +  Towsr - Spring-Dashpot 


Figure  43.  Displacements  at  the  top  of  the  tower  for  the  dam  plus  tower  and  the  frame  with 
spring-dashpot  models  (Anza) 
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Fifteen  additional  earthquake  time-histories  were  selected  and  applied  to  the 
frame  model  and  dam  plus  tower  finite  element  model.  Some  of  the  most  relevant 
characteristics  of  these  earthquakes  are  provided  in  Table  6.  This  information 
was  obtained  to  have  more  data  to  clarify  the  problem.  The  Fourier  spectrum  was 
calculated  to  obtain  the  dominant  period  and  frequency  of  the  motion.  The  dura¬ 
tion  of  the  strong  motion  was  also  calculated.  The  dominant  periods  and  frequen¬ 
cies,  total  duration,  and  duration  of  the  strong  motion  (Airong)  are  displayed  in  the 
table.  The  first  five  earthquakes  (in  bold  italics)  are  those  for  which  the  results  of 
the  frame  with  spring-dashpots  and  dam  plus  tower  finite  element  model  do  not 
match  very  well.  For  the  other  10  earthquakes  listed  in  the  table,  the  results  are  in 
good  agreement.  At  first  glance,  the  corrunon  characteristics  of  the  first  five 
earthquakes  are  their  short  durations  (less  than  1 1  sec  except  for  the  fifth  one). 
Moreover,  some  of  them  have  high  dominant  frequencies.  The  second  observa¬ 
tion  would  be  that  perhaps  the  problem  is  the  frequency  content  of  the  accelera¬ 
tion  time-history.  Therefore,  a  series  of  tests  using  seismic  input  sinusoidal 
frmctions  was  performed  to  obtain  some  patterns  and  observe  the  trends.  These 
results  are  not  presented  here  for  the  sake  of  brevity.  It  was  found  that  the 
sinusoidal  frmctions  with  the  higher  dominant  frequencies  led  to  less  similarity 
between  the  responses  calculated  with  the  two  models. 


Table  6 

Strong  Motion  Duration,  Predominant  Periods,  and  Frequencies 
for  Several  Earthquakes 

Name 

Duration 

Dstrong 

©p 

sec 

sec 

rad/sec 

Round  Vaiiey 

6.85 

3.605 

0.0811 

77.4745 

Anza 

10.32 

2.10 

0.0786 

79.9387 

Coalinga 

10.30 

2.34 

0.1326 

47.3845 

Anzap 

10.00 

1.84 

0.3567 

17.6148 

Whittier  Narrows 

20.00 

4.42 

0.2447 

25.6771 

Friuli 

22.00 

4.48 

0.4073 

15.4264 

Trinidad 

22.00 

12.17 

0.3191 

19.6903 

Imperial  Valley 

30.00 

22.86 

0.4122 

15.2431 

Pasadena 

30.00 

14.06 

0.2090 

30.0631 

El  Centro 

31.00 

23.86 

0.6035 

10.4112 

Borrego 

35.00 

25.12 

1.6614 

3.7819 

Loma  Prieta  at  Gilroy  Building 

40.00 

8.91 

1.3320 

4.7171 

Loma  Prieta  at  Anderson  Dam 

40.00 

10.51 

0.5506 

11.4111 

Cape  Mendocino 

44.00 

18.66 

2.4544 

2.5600 

Central  California 

40.00 

25.20 

0.8009 

7.8452 

Based  on  these  results,  the  first  conclusion  obtained  was  that  the  stiffness  and 
damping  coefficients  of  the  springs  and  dashpots  needed  verification.  It  is 
recalled  that  the  values  of  these  coefficients  are  only  the  function  of  the  soil 
properties  and  main  dimensions  of  the  foundation.  However,  in  theory  these 
values  also  depend  on  the  frequency  of  the  excitation,  even  though  for  simplicity 
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this  fact  is  often  neglected  in  practice.  Therefore,  it  was  decided  to  use  new  coef¬ 
ficients  that  depend  on  the  frequency  of  the  seismic  motion.  These  were  calcu¬ 
lated  using  the  expressions  developed  by  Gazetas  et  al.  (1985).  The  equations 
proposed  by  Gazetas  et  al.  (1985)  are  very  similar  to  those  used  in  the  first 
analysis,  but  they  are  multiplied  by  a  dynamic  stiffness  coefficient  that  depends 
on  the  excitation  frequency.  For  nonharmonic  excitations  such  as  those  due  to 
real  earthquakes,  it  is  necessary  to  use  the  Frequency  Response  Method  to  fully 
account  for  the  frequency  dependency  of  the  coefficients.  However,  the  variation 
with  frequency  is  not  large,  and  it  was  deemed  sufficient  to  set  the  frequency 
equal  to  the  dominant  fi'equency  of  the  ground  motion. 

The  relative  displacement  time-history  (plotted  for  only  the  Coalinga  earth¬ 
quake,  for  the  sake  of  brevity)  is  shown  as  Figure  44.  This  figure  shows  the 
responses  obtained  with  the  frame  with  spring-dashpots  in  which  the  coefficients 
are  frequency  independent  and  dependent  on  the  frequency  (the  latter  case  is 
labeled  as  new  spring-dashpot  in  the  plot).  These  displacements  are  compared 
with  those  from  the  dam  plus  tower  finite  element  model,  which  is  regarded  as 
the  “exact”  case.  This  plot  shows  no  significant  differences  in  the  results  between 
the  two  sets  of  springs  and  dashpots,  indicating  that  disregarding  the  frequency 
dependency  of  the  stiffness  and  damping  coefficients  is  not  the  cause  of  the  prob¬ 
lem.  Similar  results  were  obtained  for  the  other  earthquakes,  confirming  again 
that  there  is  no  need  to  use  frequency-dependent  stiffhess  and  dashpots  coeffi¬ 
cients,  at  least  for  this  particular  system.  This  statement  should  be  reviewed  on  a 
case-by-case  basis,  however.  One  of  the  reasons  why  there  is  no  need  to  consider 
frequency-dependent  functions  could  be  that  the  foundation  block  is  very  rigid 
and  may  not  be  sensitive  to  high-frequency  excitation  (Gazetas  1991). 


Figure  44.  Displacements  obtained  with  the  dam  plus  tower  and  spring-dashpot  models  with 

frequency  independent  and  dependent  coefficients  (Coalinga)  [Note:  -  =  Dam 

+  Tower,  —  =  Spring-Dashpot,  and - 1 - -  Spring-Dashpot-new] 
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New  sets  of  tests  were  performed  to  achieve  good  agreement  between  the 
response  obtained  with  the  simplified  and  fidi  finite  element  model  and  to  arrive 
at  conclusions  regarding  the  reasons  for  the  discrepancies. 

First,  an  eigenvalue  analysis  of  the  fixed-base  intake  tower  modeled  as  a 
frame  with  and  without  the  bridge  was  undertaken  to  observe  the  effect  of  the 
latter  on  the  frequencies  and  periods.  Next,  the  frame  with  spring-dashpots  repre¬ 
senting  a  flexible  soil  with  and  without  the  bridge  was  analyzed  to  obtain  its  nat¬ 
ural  frequencies  and  periods.  Table  7  shows  the  natural  periods  for  the  simplified 
model  of  the  intake  tower  with  and  without  the  service  bridge  and  resting  in  a 
fixed  and  flexible  soil  foundation.  Table  8  displays  the  natural  frequencies  for  the 


Table  7 

Natural  Periods  of  the  Simplified  intake  Tower  Model 


Fixed/Bridge 


Plexible/Bridge  |  Fixed/No  Bridge  |  Flexible/No  Bridge 


Table  8 

Modal  Frequencies  for  the  Simplified  Intake  Tower  Model 


Flxed/Brldge 


Frequency 


Flexible/Bridge 


Frequency 


Fixed/No  Bridge 


Frequency 


Flexibie/No  Bridge 


Frequency 

rad/sec 


49.7275 

178.7080 

268.6613 

576.1425 


1001.2780 

1419.8940 

1489.9280 

1816.4410 

2088.3500 


4.5338 

17.9409 

23.6393 

50.8007 

123.5330 

249.0922 

319.6271 

691.2986 

695.8807 

1097.0620 


48.9366 

179.4293 

268.8030 

577.0503 

755.5436 

1001.6660 

1419.9200 

1491.2250 

1816.7210 

2096.8810 


3.5268 

17.9774 

23.6270 

249.2425 

320.4843 

692.0222 

695.8858 

1097.3350 

1387.7560 

1511.9260 


same  cases.  The  cases  on  rigid  foundation  were  analyzed  to  obtain  a  measure  of 
the  flexibility  of  the  springs  on  the  spring-dashpot  or  flexible  foundation  model. 
It  can  be  noticed  that,  when  a  fixed-base  condition  is  used,  the  presence  of  the 
bridge  has  almost  no  effect  on  the  periods  and  frequencies.  For  the  flexible 
foundation  models,  the  periods  increased  as  expected  when  compared  with  the 
fixed-base  cases.  It  can  be  noticed  here  a  small  increase  in  the  lower  periods  of 
the  model  of  the  tower  without  the  bridge.  These  results  seem  to  indicate  that  the 
springs  at  the  base  of  the  tower  are  very  rigid  in  all  or  some  directions,  since  a 
greater  increase  in  the  periods  for  the  flexible-base  case  was  expected.  This  is  so 
because  the  dominant  vibration  periods  in  the  response  calculated  with  the  finite 
element  model  of  the  dam  plus  tower  (shown  in  Figures  42  and  43)  are  larger 
than  those  of  the  spring-dashpots  models.  To  verify  this  preliminary  conclusion, 
the  periods  and  frequencies  of  the  dam  plus  intake  tower  finite  element  model 
were  computed. 

Table  9  summarizes  the  natural  periods  and  frequencies  obtained  with  the 
finite  element  model  of  dam  plus  intake  tower  with  and  without  the  access 
bridge.  These  results  should  be  compared  with  the  previously  obtained  periods 
and  frequencies  corresponding  to  the  simplified  intake  tower  model  resting  on 
flexible  foundation  (the  spring-dashpot  model)  shown  in  Tables  7  and  8.  Table  9 
also  contains  natural  periods  and  frequencies  associated  with  the  vibration  of  the 
dam  itself;  thus,  these  results  cannot  be  compared  on  a  one-to-one  basis  with 
those  in  Tables  7  and  8.  Rather,  the  comparison  should  be  made  examining  the 
mode  shapes.  Nevertheless,  in  general,  the  results  confirm  that  the  finite  element 
model  consisting  of  the  dam  and  tower  is  more  flexible  than  the  spring-dashpot 


Table  9 

Natural  Periods  and  Frequencies  for  the  Dam  Plus  Tower  Model 


I  Bridge 

No  Bridge  || 

L_ 

i 

mim 

1 

1 

1  111= 

i  1  1  F 

Mode 

No  Bridge 

Bridge 

Mode 

No  Bridge 

Bridge 

Period 

sec 

Period 

sec 

Frequency 

rad/sec 

Frequency 

rad/sec 

1 

19.7359 

2.3753 

1 

0.3184 

2.6452 

2 

1.6888 

1.7573 

2 

3.7204 

3.5755 

3 

1.2536 

1.2940 

3 

5.0122 

4.8556 

4 

0.6309 

0.5776 

4 

9.9588 

10.8776 

5 

0.5643 

0.5642 

5 

11.1349 

11.1362 

6 

0.5445 

0.5432 

6 

11.5399 

11.5674 

7 

0.5165 

0.5155 

7 

12.1647 

12.1896 

8 

0.5050 

0.5055 

8 

12.4424 

12.4304 

9 

0.4566 

0.4566 

9 

13.7603 

13.7611 

10 

0.4312 

0.4340 

10 

14.5724 

14.4764 
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model.  Note  also  that  when  the  bridge  is  taken  out,  the  periods  increase  more 
dramatically  than  in  the  case  of  the  spring-dashpot  model.  These  findings  demon¬ 
strate  that,  for  the  present  system,  the  values  of  the  spring  constants  used  to  rep¬ 
resent  the  flexible  soil  foundation  are  not  adequate  to  get  a  good  approximation 
with  the  simple  model.  They  underestimate  the  flexibility  of  the  soil. 

Before  proceeding  to  deal  with  this  situation,  it  is  very  instructive  to  examine 
the  mode  of  vibration  corresponding  to  the  lowest  natural  period  of  the  two  finite 
element  models  considered  in  Table  9.  These  are  shown  in  Figures  45  and  46. 
Note  that,  as  shown  in  Figure  46,  the  very  high  period  (19.7  sec)  corresponds  to  a 
rigid  body  rotation  of  the  intake  tower  caused  by  the  rocking  flexibility  of  the 
soil.  The  situation  is  similar  when  the  bridge  is  included  in  the  model  (Figure  45), 
but  the  value  of  the  period  is  now  restricted  by  the  presence  of  the  bridge. 


Figure  45.  Mode  shape  obtained  with  the  dam  plus  tower  and  bridge  finite  element  model 


Figure  46.  Mode  shape  obtained  with  the  dam  plus  tower  finite  element  model  without  the 
bridge 


It  was  decided  to  evaluate  the  stiffiiess  coefficients  of  the  springs  at  the  tower 
base  using  another  method.  Ideally,  the  coefficients  should  be  determined  using 
an  analytical  method  so  that  they  can  be  later  applied  to  other  situations. 
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However,  this  is  a  task  well  beyond  the  scope  of  this  report.  Moreover,  it  is  not 
clear  what  assumptions  used  to  define  the  coefficients  currently  available  are  not 
appropriate  and  cause  the  problem  previously  described.  Therefore,  it  was 
decided  to  perform  a  numerical  simulation.  As  was  done  to  determine  the  springs 
at  the  support  of  the  access  bridge,  the  stiffness  coefficients  were  calculated 
using  the  unit  force  method  and  SAP2000.  Table  10  shows  the  values  for  the 
stiffness  coefficients  obtained  with  this  approach.  The  unit  load  was  applied  in 
the  horizontal  and  vertical  directions,  and  a  unit  moment  was  applied  along  the 
axis  perpendicular  to  the  plane  of  analysis.  The  values  identified  as  “Equations” 
correspond  to  those  obtained  using  the  formulas  in  Chapter  4;  i.e.,  they  were 
calculated  with  the  equations  published  in  the  literature.  These  values  were  used 
in  the  first  spring-dashpot  model.  and  Ky  represent  the  stiffness  of  the  soil  in 
the  horizontal  and  vertical  directions,  respectively,  and  Kr  is  the  coefficient  of  the 
rocking  spring.  It  can  be  noticed  that  the  coefficients  calculated  with  the  two 
methods  show  substantial  differences. 


Table  10 

Stiffness  Coefficients  for  the  Spring-Dashpot  Model 

K, 

Ky 

Kr 

Method 

kip/ft 

kip/ft 

kip-ft/rad 

Equations 

47,367 

49,691 

6,245,753 

Unit-Force 

1,265 

2,406 

44,763 

Several  tests  were  done  varying  the  three  original  stiffiiess  coefficients.  The 
results  demonstrated  that  the  major  factor  influencing  the  response  of  this  struc¬ 
tural  system  is  the  rocking  stiffness.  The  best  agreement  between  the  dam  plus 
tower  finite  element  model  and  the  frame  with  spring-dashpots  was  obtained 
when  the  rocking  stiffness  calculated  from  the  unit-force  method  was  used.  The 
periods  and  frequencies  from  the  eigenvalue  analysis  were  calculated  and  exam¬ 
ined  for  a  second  time.  The  first  10  new  natural  periods  and  frequencies  are  pre¬ 
sented  in  Tables  1 1  and  12,  respectively.  These  should  be  compared  with  those  in 
Table  9.  It  is  evident  that,  using  the  new  stiffness  coefficients,  the  periods  and 
fi'equencies  of  the  frame  with  spring-dashpots  are  much  closer  to  those  of  the 
dam  plus  tower  finite  element  model.  The  natural  periods  obtained  with  the  new 
rocking  coefficient  in  Table  10  increased  and  are  now  closer  to  the  corresponding 
periods  of  the  dam  plus  tower  model. 

Once  the  natural  periods  were  in  better  agreement,  the  displacements  of  the 
intake  tower  with  the  new  stiffness  coefficients  were  calculated  again.  The  results 
of  the  new  time-history  analyses  are  presented  in  Figures  47  through  52.  These 
figures  show  a  comparison  between  the  relative  displacements  at  the  top  of  the 
tower  obtained  by  using  the  dam  plus  tower  finite  element  model  and  the  frame 
with  spring-dashpots.  The  tower  response  is  presented  for  several  earthquakes 
with  different  durations  and  frequency  contents.  For  a  same  earthquake,  the 
displacements  were  calculated  in  two  ways:  first,  all  the  vibration  modes  were 
used  to  calculate  the  response;  next,  the  response  was  obtained  by  considering 
only  the  contribution  of  the  first  mode. 
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Figure  47.  Relative  dispiacements  obtained  with  the  dam  pius  tower  and  the  new  spring- 
dashpot  models  including  all  modes  (Anza) 


Figure  48.  Relative  dispiacements  obtained  with  the  dam  plus  tower  and  the  new  spring- 
dashpot  models  considering  only  the  first  mode  of  vibration  (Anza) 
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Displacement  (in.)  ^  Displacement  (in.) 


3  49.  Relative  displacements  obtained  with  the  dam  plus  tower  and  the  new  spring- 
dashpot  models  including  all  modes  (Coalinga) 


Figure  50.  Relative  displacements  obtained  with  the  dam  plus  tower  and  the  new  spring- 
dashpot  models  considering  only  the  first  mode  of  vibration  (Coalinga) 
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re  51 .  Relative  displacements  obtained  with  the  dam  plus  tower  and  the  new  spring 
dashpot  models  including  all  vibration  modes  (El  Centro) 


Figure  52.  Relative  displacements  obtained  with  the  dam  plus  tower  and  the  new  spring 
dashpot  models  considering  only  the  first  mode  of  vibration  (El  Centro) 
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The  first  four  plots  (Figures  47  though  50)  correspond  to  the  short-duration 
earthquakes  (Anza  and  Coalinga).  These  graphs  are  presented  first,  since  it  was 
for  these  earthquakes  that  the  larger  differences  between  the  results  of  the  spring- 
dashpot  model  and  the  finite  element  analysis  were  observed.  The  last  two  plots 
(Figures  51  and  52)  show  the  displacements  at  the  top  of  the  tower  for  the  El 
Centro  earthquake,  for  which  reasonable  results  were  obtained  with  the  original 
spring  coefficients.  A  better  agreement  between  the  two  models  was  now 
achieved. 

It  is  evident  that  now  all  the  plots  show  a  good  agreement  between  the 
responses  calculated  with  both  models.  This  fact  verifies  the  importance  of  the 
soil  rocking  stiffiiess  in  the  seismic  analysis  of  intake  towers  founded  in  soft 
soils.  Moreover,  the  results  indicate  the  need  of  refining  the  definitions  of  the 
spring  coefficients  so  that  they  more  accurately  represent  the  stiffiiess  of  the  soil. 

Comparing  the  time-histories  computed  with  several  modes  and  with  the 
lowest  mode,  one  can  conclude  that  a  good  approximation  can  be  achieved  by 
using  only  the  first  mode  of  vibration.  This  occurs  because  the  soil-structure 
interaction  accentuates  the  importance  of  the  first  mode  of  vibration,  i.e.,  the 
rigid  body  rotation  of  the  tower. 

The  last  set  of  results  demonstrates  that  the  proposed  fi^ame  with  spring- 
dashpots  model  is  acceptable  for  calculating  the  approximate  seismic  response  of 
the  intake  tower  and  bridge  system,  provided  that  the  proper  spring  coefficients 
are  adopted.  Note  that  the  frame  model  slightly  overestimates  the  magnitude  of 
the  displacements.  These  differences  could  be  lessened  by  changing  the  values  of 
the  damping  coefficients.  However,  in  this  study,  the  coefficients  of  the  dashpots 
were  not  modified;  i.e.,  the  original  expressions  in  Table  5  were  used  to  define 
them. 

The  importance  of  the  rocking  stiffness  in  this  case  study  could  be  the  result 
of  many  factors.  Some  studies  confirm  that,  for  tall  structures  such  as  the  pylons 
of  tall  bridges,  the  most  important  mode  of  foundation  behavior  is  the  rocking 
stiffness  term  (Po  Lam  and  Law  2000).  The  intake  tower  of  this  study  is  a 
relatively  tall  structure  founded  on  soft  soil  causing  the  increase  of  the  displace¬ 
ments  and  rotation  at  its  base.  Thus,  perhaps  for  this  reason,  the  rocking  stiffness 
has  a  great  influence  in  the  overall  response  of  this  particular  case  study.  Also, 
with  regard  to  the  subject  of  rocking  in  buildings,  another  author  (Bard  1988) 
conducted  a  study  measuring  the  transverse  and  rocking  components  of  the 
motion  at  the  base  of  several  tall  buildings.  The  buildings  of  this  study  have 
between  3  and  12  stories.  Bard  (1988)  confirmed  that  if  the  building  is  founded 
on  a  flexible  base,  the  rocking  motion  should  be  larger  than  for  the  same  building 
founded  on  a  fixed  base.  He  stated  that  the  large  amount  of  rocking  proved  the 
existence  of  a  very  strong  soil-structure  interaction.  These  two  studies  verified 
the  importance  of  having  a  good  approximation  of  the  rocking  stiffness  to  prop¬ 
erly  account  for  the  soil-structure  interaction  phenomenon. 


78 


Chapter  6  Numerical  Results 


Response  Spectrum  Analysis 

A  response  spectrum  analysis  was  performed  to  directly  obtain  the  maximum 
responses  of  the  intake-outlet  tower  under  study.  The  base  shear  and  overturning 
moment  were  calculated  by  implementing  the  formulation  described  below  into  a 
MATLAB-based  program. 

First,  the  maximmn  modal  forces  for  a  frame  model  were  calculated  as 
follows: 


T 

--K0,.  " 

max  J 


,  PSAj 


(45) 


where 

K  =  stiffness  matrix 

0j  =  mode  of  vibration  column  of  the  eigenvector  matrix) 

Tj  =  corresponding  period 

PSAj  =  spectral  acceleration  obtained  from  the  response  spectrum  for  the 
period  Tj 

For  a  plane  frame  model.  Equation  45  gives  forces  in  the  horizontal  and 
vertical  directions  as  well  as  moments.  However,  if  a  lumped  mass  model  is  used, 
the  rotations  (and  moments)  are  condensed  out;  i.e.,  only  the  translational  dofs 
are  used.  Adding  the  horizontal  components  of  the  forces  in  Equation  45,  one 
obtains  the  modal  base  shear  for  a  given  mode  j  with 


■j  =  \..N 


(46) 


in  which  N  is  the  total  number  of  dofs.  In  Equation  6.3  it  is  assumed  that  the 
horizontal  forces  are  those  with  odd  subindices,  and  a  lumped  mass  model  was 
used. 

The  total  maximum  base  shear  can  be  approximately  computed  by  using  the 
SRSS  (square  root  of  the  sum  of  the  squares)  modal  combination  rule: 


v  =  ^  (47) 

The  modal  overturning  moments  are  calculated  using  the  horizontal  modal 
forces: 


N 


Fh  ;;=1...A^ 

ij  i 


(48) 
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where  hi  is  the  distance  from  the  point  of  application  of  the  force  Fjj  to  the  base. 
Using  the  SRSS  modal  combination  rule,  the  total  maximiun  overturning 
moment  at  the  base  is  obtained  as  follows: 


m=tm; 

Vi=i 


(49) 


The  base  shear  and  overturning  moment  were  calculated  using  a  design 
spectrum  for  a  soil  type  Sj)  and  damping  ratio  of  24,5  percent.  It  is  recalled  that 
this  is  the  equivalent  damping  ratio  that  accounts  for  the  nonlinear  behavior  of 
the  soil  and  was  obtained  by  matching  the  responses  from  QUAD4M  and 
SAP2000.  These  results  were  compared  with  those  obtained  from  a  time-history 
analysis  using  earthquakes  compatible  with  the  design  spectrum.  To  obtain  a 
design  spectrum  for  24.5  percent  damping  ratio,  the  reduction  factors  developed 
by  Newmark  and  Hall  (1982)  were  applied  to  the  spectral  acceleration  spectrum 
for  5  percent  of  damping.  The  reduction  factor  used  is 


SR 

A 


3.21-0.68ln(f^^) 

2A2 


(50) 


where  Ceji^  the  desired  damping  ratio  (in  this  case,  0.245)  and  the  constant 
“2.12’^  in  the  denominator  is  the  corresponding  amplification  factor  with  respect 
to  the  5  percent  damped  elastic  design  spectrum.  The  design  spectrum  selected 
for  the  study  was  the  one  specified  in  the  1997  Uniform  Building  Code  (ICBO 
1997)  for  seismic  zone  3  and  soil  profile  type  So. 

A  set  of  compatible  earthquakes  was  produced  by  Luis  A.  Montejo,  graduate 
student  at  the  University  of  Puerto  Rico  at  Mayagiiez.  He  generated  two  compati¬ 
ble  earthquakes  for  a  UBC-97  design  spectrum  with  24.5  percent  damping  and 
soil  type  The  process  for  the  generation  of  compatible  accelerograms  is  based 
on  the  modification  of  recorded  acceleration  time-histories  using  the  wavelet 
transform  (Montejo  2004).  The  El  Centro  and  Loma  Prieta  earthqxiakes  were 
used  as  starting  or  seed  accelerograms  to  implement  the  process.  In  addition  to 
the  analysis  using  the  UBC-97  design  spectrum  and  to  the  spectrum-compatible 
accelerograms,  a  response  spectrum  analysis  using  the  ground  response  spectra 
of  the  El  Centro  and  the  Loma  Prieta  earthquakes  was  undertaken. 

Figures  53  and  54  show  the  pseudo-acceleration  response  spectra  for  the 
El  Centro  and  Loma  Prieta  earthquakes,  respectively.  Each  figure  contains  three 
spectra:  the  5  percent  damping  spectrum  of  the  original  earthquake,  the  24.5  per¬ 
cent  of  the  original  accelerogram,  and  die  spectrum  compatible  with  the  UBC-97 
for  24.5  percent  damping. 

Tables  13  through  16  present  the  base  shear  and  overturning  moment  calcu¬ 
lated  with  the  response  spectrum  or  time-history  analysis.  Two  models  of  the 
tower  were  used:  the  frame  model  previously  described  and  a  shear  building 
model.  Shear  building  models  are  not  commonly  used  to  study  structures  such  as 
towers  because  the  deformation  pattern  of  the  latter  structural  systems  is  gov¬ 
erned  by  bending.  However,  it  was  considered  interesting  to  examine  how  well 
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Figure  53.  Acceleration  response  spectra  for  the  El  Centro  earthquake 
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Figure  54.  Acceleration  response  spectra  for  the  Loma  Prieta  earthquake 


Table  13  I 

Base  Shears  and  Overturning  Moments  for  El  Centro  and  Loma 
Prieta  Earthquakes  (Response  Spectrum  Analysis)  I 

Earthquake 

Type  of  Analysis 

Shear  Building  Model 

Frame  Model  | 

Vi,ase 

kip 

Mbase 

kip-ft 

V,ase 

kip 

Mbase 

kip-ft 

Loma  Prieta 

Spectrum  (24.5  %) 

1445.50 

35315.00 

1442.90 

24675.00 

El  Centro 

Spectrum  (24.5  %) 

1300.00 

28159.00 

1597.00 

26664.00 

I  Table  14 

Base  Shears  and  0\ 
1  Prieta  Earthquakes 

rerturning  Moments  for  El  Centro  and  Loma 
[Time-History  Analysis)  I 

1  Earthquake 

Type  of  Analysis 

Frame  Model  1 

\4ase 

kip 

Afftase 

kip-ft 

1  Loma  Prieta 

Time-history 

1099.25 

39379.25 

1  El  Centro 

Time-history 

1222.78 

30033.13 
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Table  15 

Base  Shears  and  Overturning  Moments  for  El  Centro  and  Loma 
Prieta  Spectrum-Compatible  Accelerograms  (Response  Spectrum 
Analysis) 

Earthquake 

Type  of  Analysis 

Shear  Building  Model 

Frame  Model  | 

Vi.„ 

kip 

Mtase 

kip-ft 

14a,. 

kip 

Mbase 

klp.ft 

Loma  Prieta 

Spectrum  (24.5  %) 

1486.40 

35911.00 

1732.20 

28845.00 

E(  Centro 

Spectrum  (24.5  %) 

1474.00 

34297.00 

1849.80 

29984.00 

Table  16 

Base  Shears  and  Overturning  Moments  for  El  Centro  and  Loma 
Prieta  Spectrum-Compatibie  Accelerograms  (Time-History 

Anaiysis) 

Earthquake 

Type  of  Analysis 

Frame  Model  I 

Vt». 

kip 

Mtase 

kip-ft 

Loma  Prieta 

Time-history 

1293.13 

39862.25 

El  Centro 

Time-history 

1548.90 

29426.00 

the  shear  building  model  approximated  the  seismic  response  of  intake  towers. 
The  advantage  of  this  model  is  its  simplicity. 

The  agreement  between  the  results  of  the  time-history  and  spectrum  analyses 
is  quite  reasonable  for  both  earthquakes,  and  for  the  original  and  spectrum- 
compatible  accelerograms.  It  is  interesting  to  observe  from  the  values  in  the 
tables  that  the  agreement  between  the  time-history  and  spectrum  analyses  is  bet¬ 
ter  when  the  shear  building  model  is  used.  This  could  be  because  the  errors  intro¬ 
duced  when  a  shear  building  model  is  used  are  somehow  compensating  for  the 
approximations  inherent  in  the  response  spectrum  approach. 


Summary 

This  chapter  presented  numerical  examples  of  the  seismic  analysis  of  intake- 
outlet  towers  and  their  evaluation.  The  numerical  tools  used  for  the  analysis  of 
the  towers  that  accounts  for  the  multiple-support  excitation  and  soil-structure 
interaction  effects  were  two  finite  element  programs  and  a  MATLAB-based 
program.  A  typical  intake-outlet  tower  from  the  USAGE  inventory  was  selected 
for  all  the  analyses.  Thus,  this  study  can  help  to  establish  a  behavioral  pattern  of 
the  dynamic  response  of  comparable  towers  founded  on  similar  soil  conditions. 

The  response  of  the  towers,  including  the  access  bridge  and  multiple-support 
excitations,  indicates  that  these  effects  are  not  significant,  at  least  for  this  case 
study.  The  amplification  of  the  ground  acceleration  at  the  bridge  support  was  not 
significant,  probably  due  to  the  particular  soil  conditions  and  geometry  of  the 
dam.  Another  important  finding  of  the  study  is  that  the  dynamic  response  of  this 
intake  tower  is  mainly  influenced  by  the  soil-structure  interaction  effects.  The 
soil-structure  interaction  in  this  system  was  very  strong  due  to  the  large 
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differences  between  the  stiffhess  of  the  soil  and  the  tower.  That  is,  the  intake 
tower  is  very  rigid  compared  with  the  soft  soil  beneath  it.  The  results  also 
indicate  that  the  rocking  stif&iess  is  the  parameter  that  plays  the  most  important 
role  in  the  dynamic  response  of  the  structural  system  studied.  Therefore,  it  is 
crucial  to  have  accurate  expressions  to  evaluate  the  stiffiiess  of  the  spring 
associated  with  the  rocking  motion. 

The  values  obtained  with  commonly  used  formulas  did  not  provide  satisfac¬ 
tory  results  in  the  present  study.  To  overcome  this  problem,  new  values  of  the 
stifihess  coefficients  were  obtained  by  applying  a  unit  force  (or  moment)  to  the 
finite  element  model  at  the  location  of  the  foundation  and  calculating  the  dis¬ 
placement  (or  rotation).  It  is  realized,  however,  that  this  approach  is  not  practical 
and  thus  more  research  is  needed  on  this  subject.  Finally,  it  is  important  to 
emphasize  ftiat  the  effects  of  multiple-support  excitation  and  soil-structure 
interaction  should  be  reviewed  on  a  case-by-case  basis.  It  may  happen  that,  for 
other  towers  and  dams,  the  importance  of  these  two  effects,  or  the  lack  of  it, 
could  be  reversed. 
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7  Conclusions  and 
Recommendations 


Summary  and  Conclusions 


Hie  intake-outlet  was  selected  based  on  data  available  from  the  USAGE  mvent- 
toiy.  The  Aikabutla  Dam  located  in  northern  Mississippi,  30  miles  south  of 
Memphis,  TN,  was  chosen  for  the  analysis  because  it  is  of  particular  interest  for 
the  Corps  of  Engineers.  Hie  geometry  of  the  original  intake  tower  was  changed 
so  that  it  became  relatively  simple  to  permit  a  simplification  of  the  analysis.  TTie 
analj^ses  were  performed  using  severd  two-dimensional  finite  element  models 
and  a  fiame  model.  A  computer  program  w^  also  developed  for  this  investiga¬ 
tion  using  MATLAB.  Two  finite  element  computer  programs  were  used: 
QUAD4Mand  SAP2000.  The  finite  element  models  developed  are  the  Aikabutla 
Dam  (using  QUAD4Mand  SAP2000)  and  the  Arkabutla  Dam  plus  fire  intake 
tower  (using  SAP2000).  The  models  implemented  in  the  MATLAB  program  are 
tv^'O  plane  flume  models  of  the  intake  tower  with  and  without  the  bridge,  and  the 
same  systems  widi  springs  and  dashpots  at  the  supports.  Included  in  the 
MATLAB  program  are  the  formulation  for  multiple-support  excitetion  effects, 
the  soil-structure  interaction  problem,  and  the  simplified  procedure  for  the  calcu- 
at 


models.  The  main  response  quantity  used  for  the  companson  is  the  relative  dis¬ 
placement  at  the  top  of  the  tower.  The  results  are  presented  in  terms  of  plots  of 
the  top  relative  displacement  as  fimction  of  time.  The  base  shear  and  overturning 
moments  were  calculated  when  a  response  spectrum  analysis  was  done. 

One  of  the  principal  findings  of  this  report  was  that  the  effect  of  the  multiple- 
support  excitation  in  the  response  of  the  tower  and  service  bridge  is  not  important 
and  it  can  be  disregarded.  Ihe  slight  variation  in  the  soil  properties  between  the 
base  of  the  tower  and  the  location  of  the  bridge,  combined  with  low  ^pect  ratio 
of  earth  dam  eeometr%'  may  have  contributed  to  this  result.  Due  to  the  similar 


Chapter  7  Conclusions  ar 


expected  that  there  would  be  significant  differences  between  the  acceleration 
time-histories  at  the  two  supports,  but  the  results  demonstrated  otherwise. 

Another  important  result  obtained  from  this  research  is  the  demonstration  of 
the  important  influence  of  the  soil-structure  interaction  effects  in  the  response  of 
the  tower.  Initially,  this  effect  was  not  going  to  be  included  in  the  study  because 
it  was  not  deemed  to  be  important.  However,  the  completely  different  responses 
obtained  with  the  model  of  the  tower  plus  bridge  with  fixed  bases  and  those  from 
the  “exact”  finite  element  model  led  to  the  authors  investigating  this  effect.  By 
examining  the  response  of  the  tower  from  the  finite  element  model,  it  could  be 
seen  that  the  displacements  were  the  sum  of  those  due  to  the  bending  deforma¬ 
tion  plus  those  due  to  the  rotation  of  the  tower  as  a  rigid  body.  To  confirm  that 
the  soil-structure  interaction  was  responsible  for  the  differences  observed,  the 
accelerations  at  the  interface  between  the  tower  and  the  soil  with  the  tower 
included  in  the  finite  element  model  were  compared  with  those  obtained  at  the 
free  surface  of  the  model  of  the  dam  without  the  tower.  As  was  suspected,  the 
accelerations  were  quite  different. 

It  was  then  decided  to  accormt  for  the  soil-stmcture  interactions  by  simple 
means  so  that  they  can  be  incorporated  in  the  frame  model  of  the  tower  and 
bridge  system.  The  simplest  way  to  consider  the  soil-stmcture  interaction  is  by 
means  of  springs  and  dashpots  at  the  base  of  the  stmcture.  This  idea  was  origi¬ 
nally  developed  for  the  study  of  vibration  of  machine  foundations.  Then  the 
following  problem  emerged:  the  formulas  that  define  the  stiffiiess  and  damping 
coefficients  of  the  springs  and  dashpots  carmot  be  applied  at  the  end  of  the  bridge 
that  lies  on  the  dam.  Indeed,  these  coefficients  were  obtained  under  the  premise 
that  a  harmonic  force  was  ^plied  to  a  rigid  circular  disk  resting  on  top  of  an 
elastic  half-space.  Clearly,  Ae  crest  of  a  dam  is  far  from  this  assumption.  There¬ 
fore,  it  was  decided  to  obtain  the  stiffiiess  coefficients  by  ^plying  a  imit  force  at 
the  top  of  the  dam,  one  at  a  time  in  two  directions  (horizontal  and  vertical). 
Moreover,  it  was  also  decided  not  to  place  any  dashpots  there  (a  conservative 
assumption). 

To  obtain  more  accurate  values  for  the  stiffiiess  coefficients,  as  well  as  to 
determine  the  dashpot  coefficients,  one  has  to  apply  a  force  of  the  form  and 

obtain  the  steady-state  response  for  different  excitation  frequencies  co.  This 
option  was  not  pursued  because  it  was  beyond  the  scope  of  the  report.  The  coef¬ 
ficients  for  the  springs  and  dashpots  in  three  directions  (horizontal,  vertical,  and 
rotation  about  an  axis  normal  to  the  plane)  at  the  tower  base  were  taken  from  the 
technical  literature.  When  the  response  of  the  frame  with  the  springs  and  dash- 
pots  was  compared  with  the  finite  element  model,  the  matching  was  acceptable 
for  certain  earthquakes  but  inaccurate  for  others.  Additionally,  a  comparison  of 
the  tower’s  natural  periods  detected  that  the  periods  of  the  finite  element  model 
were  larger  than  those  predicted  by  the  frame  with  springs  and  dashpots.  By 
varying  the  values  of  the  spring  coefficients  to  try  to  obtain  a  better  agreement,  it 
was  found  that  the  coefficient  of  the  rocking  spring  was  responsible  for  the  mis¬ 
match.  This  is  another  finding  of  the  research:  at  least  for  the  case  studied  in  this 
investigation,  the  coefficient  associated  with  the  rocking  motion  of  the  founda¬ 
tion  does  not  well  represent  the  flexibility  of  the  soil.  Therefore,  it  was  concluded 
that  the  value  of  the  rocking  stiffiiess  would  have  to  be  determined  by  applying  a 
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unit  moment  at  the  free  surface.  When  used  with  die  frame  model,  the  value 
obtained  in  this  way  allowed  a  very  good  agreement  with  the  results  from  the  full 
finite  element  model  of  the  dam-tower-bridge-foundation  system. 

It  can  be  argued  that  the  soil  underneath  die  reservoir  of  the  case  studied  w'as 
soft  and  that  in  other  cases  the  soil-structure  interaction  will  not  be  important. 
However,  to  establish  if  this  is  indeed  the  case,  a  brief  parametric  study  was 
undertaken  in  which  die  stiffiiess  of  the  soil  deposit  (the  shear  wave  velocity  or 
elastic  moduli)  was  gradually  increased.  It  was  found  that  the  effect  remained 
important  except  for  a  material  that  could  be  characterized  as  a  hard  rock. 


Recommendations  for 
Further  Studies 

Although  this  study  tried  to  be  as  comprehensive  as  possible  within  the  time 
constraints,  diere  are  a  number  of  topics  or  ideas  that  were  not  addressed  or 
explored  because  of  time  limitations  or  because  they  fall  outside  the  scope  of  die 
study.  How^ever,  it  is  beheved  that  studying  these  topics  is  w'orthwhile  and  can 
complement  the  results  of  this  work.  The  subjects  recommended  for  fiirther 
studies  are  listed  below  in  no  specific  order. 

a.  A  nonlinear  analysis  of  the  intake-outlet  should  be  performed  to  examine 
the  dynamic  behavior  of  this  kind  of  stracture  until  feilure  or  when  the 
tower  begins  to  degrade  during  a  severe  earthquake. 

b.  The  nonlinear  analysis  can  be  extended  to  include  the  model  of  the  dam 
plus  the  intake  tower.  When  the  nonlinear  behavior  of  the  soil  is  consid¬ 
ered,  the  stiffiiess  of  the  soil  changes  with  time,  depending  on  die  level 
of  strains  (Kumar  and  Prakash  1997).  Thus,  at  any  particular  time,  the 
strains  in  soil  and  dius  the  foundation  stiffiiess  will  not  be  the  same  for 
any  two  different  base  input  motions.  It  will  be  very  interesting  to  study 
this  phenomenon  for  intAe-outlet  towers.  To  do  that,  it  will  be  necessary 
to  use  another  computer  program  with  the  capabilities  to  perform 
nonlinear  analysis  of  soil  stractures  in  addition  to  conventional  ones. 

c.  The  same  methodology  used  in  this  study  can  be  applied  to  other  towers 
with  typical  dimensions  but  for  taller  dams.  The  multiple-support  excita¬ 
tion  effects  could  be  important  in  these  cases,  especially  if  the  material  of 
the  dam  is  quite  different  from  the  soil  beneath  the  tower. 

d  The  accuracy  of  the  stiffness  coefficients  of  the  springs  representing  die 
rocking  motion  should  be  examined.  To  properly  account  for  the  soil- 
stracture  interaction  effects  in  tall,  slender  structures  such  as  intake 
towers  and  chimneys,  it  is  cracial  to  have  rchable  values  of  these 
coefficients. 

e.  It  will  be  useful  for  future  analysis  of  tower-bridge  systems  to  have 
spring  and  dashpot  coefficients  obtained  for  a  wedge-like  soil  structure 
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such  as  an  earthen  dam.  This  could,  in  principle,  be  done  in  an  anal5l;ical 
fashion  following  a  procedure  similar  to  that  for  the  elastic  half-space 
case.  However,  due  to  the  complexity  of  the  geometry,  it  is  more  likely 
that  the  study  will  have  to  be  done  numerically.  In  this  case,  many  geo¬ 
metrical  configurations  of  the  dam  will  have  to  be  studied  in  order  to 
obtain  coefficients  that  are  a  fimction  of  the  aspect  ratio  of  the  dam 
(height/base  length)  and  the  properties  of  the  predominant  soil  type  (G, 
p,  and  u). 

f  The  survey  conducted  by  the  U.S.  Army  Corps  of  Engineers  in  which 
they  identified  the  typical  properties  of  the  intake-outlet  towers  under 
their  jurisdiction  should  be  extended  to  include  the  typical  configurations 
of  the  service  bridges.  An  inventory  of  the  service  bridges  could  be  very 
usefiil  in  identifying  the  typical  bridge  supports,  sections,  and  other 
details  that  are  necessary  in  the  modeling  process  of  this  kind  of  system. 

g.  It  would  be  interesting  to  include  the  effects  of  the  conduits  or  tunnels  in 
the  seismic  analysis  of  the  intake-outlet  towers.  These  conduits  or  tun¬ 
nels  are  very  large  and  could  play  an  important  role  in  the  response  of  the 
tower  with  the  dam. 

h.  Extending  the  analysis  of  the  tower  and  bridge  systems  in  the  perpendic¬ 
ular  direction,  i.e.,  along  the  longitudinal  axis  of  the  dam,  should  be 
done.  This  could  be  very  useful  in  identifying  the  effect  of  the  torsion 
produced  by  the  bridge  in  the  tower  and  its  influence  in  the  overall 
response  of  the  system. 

i.  Special  attention  should  be  paid  to  the  support  of  the  bridge.  Depending 
on  the  type  of  foundation,  the  bridge  might  slide  along  the  foundation  as 
the  result  of  a  strong  earthquake. 
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